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Introduction 



The Standard Model (SM) is the theory that provides the best description of the 
properties and interactions of elementary particles. The strong interaction be- 
tween quarks and gluons is described by the Quantum Chromodynamics (QCD) 
field theory. Jet production is the high-p^ process with the largest cross sec- 
tion at hadron colliders. The jet cross section measurement is a fundamental 
test of QCD and it is sensitive to the presence of new physics. It also provides 
information on the parton distribution functions and the strong coupling. One 
of the fundamental elements of jet measurements is the proper understanding of 
the energy flow around the jet core and the validation of the QCD description 
contained in the event generators, such as parton shower cascades, and the frag- 
mentation and underlying event models. Jet shapes observables are sensitive to 
these phenomena and thus very adequate to this purpose. The first measure- 
ment of the inclusive jet cross section in pp collisions at y/s = 7 TeV delivered 
by the LHC was done using an integrated luminosity of 17 nb~^ recorded by the 
ATLAS experiment. The measurement was performed for jets with Pt > 60 GeV 
and \y\ < 2.8, reconstructed with the anti-Zcj algorithm with radius parameters 
R = 0.4 and R = 0.6. 

This Ph.D. Thesis presents the updated measurement of the inclusive jet 
cross section using the full 2010 data set, corresponding to 37 pb~^ collected by 
ATLAS. Jets with pt > 20 GeV and \y\ < 4.4 are considered in this analysis. 
The measurement of the jet shapes using the first 3 pb~^ is also presented, for 
jets with Pt > 30 GeV and \y\ < 2.8. Both measurements are unfolded back to 
the particle level. The inclusive jet cross section measurement is compared to 
NLO predictions corrected for non-perturbativc effects, and to predictions from 
an event generator that includes NLO matrix elements. Jet shapes measure- 
ments are compared to the predictions from several LO matrix elements event 
generators. 

The contents of this Thesis are organized as follows: Chapter 1 contains a 
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description of the strong interaction theory and jet phenomenology. The LHC 
coUider and the ATLAS experiment are described in Chapter 2. The inclusive 
jet cross section measurement is described in detail in Chapter 3, and the jet 
shapes measurements in Chapter 4. Additional comparison of the jet shapes 
measurement to Monte Carlo event generator predictions are shown in Chapter 5. 
There are two appendixes at the end of the document. The first one contains 
additional jet shapes studies, and the second one is devoted to energy flow studies 
at calorimeter level. 
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Chapter 1 

QCD at Hadron Colliders 



1.1 The Standard Model 

The Standard Model (SM) |py is the most successful theory describing the prop- 
erties and interactions (electromagnetic, weak and strong) of the elementary 
particles. The SM is a gauge quantum field theory based in the symmetry 
group SU{3)c X SU{2)l x f/(l)y, where the electroweak sector is based in the 
SU{2)l X U{1)y group, and the strong sector is based in the SU{?>)c group. 

Interactions in the SM occur via the exchange of integer spin bosons. The 
mediators of the electromagnetic and strong interactions, the photon and eight 
gluons respectively, are massless. The weak force acts via the exchange of three 
massive bosons, the and the Z. 

The other elementary particles in the SM are half-integer spin fermions: six 
quarks and six leptons. Both interact electroweakly, but only quarks feel the 
strong interaction. Electrons (e), muons(/x) and taus(r) are massive leptons and 
have electrical charge Q = -1. Their associated neutrinos (z/g, z/^, Vr) do not have 
electrical charge. Quarks can be classified in up-type (m, s, t) and down-type 
{d,s,b) depending on their electrical charge (Q = 2/3 and Q = -1/3 respectively). 
For each particle in the SM, there is an anti-particle with opposite quantum 
numbers. 

The SM formalism is written for massless particles and the Higgs mechanism 
of spontaneous symmetry breaking is proposed for generating non-zero boson 
and fermion masses. The symmetry breaking requires the introduction of a new 
field that leads to the existence of a new massive boson, the Higgs boson, that 
has still not been observed. 
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1.2 Quantum Chromo dynamics Theory 

Quantum Chromodynamics (QCD) [2] is the renormalizable gauge field theory 
that describes the strong interaction between colored particles in the SM. It is 
based in the SU{3) symmetric group, and its lagrangian reads: 

Cqcd = -\f^pF2^ + Y: qi^YD,-m)q (1.1) 

flavors 

where the sum runs over the six different types of quarks, q, that have mass m. 
The field strength tensor, is derived from the gluon field A^'- 

F^, = [d^A^ - d,At - gf^^'^AlA^,] (1.2) 

f^^'^ are the structure constants of SU{?>)^ and the indices A, B, C run over 
the eight color degrees of freedom of the gluon field. The third term originates 
from the non-abelian character of the SU{?>) group, and is the responsible of the 
gluon self-interaction, giving rise to triple and quadruple gluon vertexes. This 
leads to a strong coupling, = 0^/4,71 that is large at low energies and small at 
high energies (see Figure fTTT]) . Two consequences follow from this: 

• Confinement: The color field potential increases linearly with the dis- 
tance, and therefore quarks and gluons can never be observed as free par- 
ticles. They are always inside hadrons, either mesons (quark-antiquark) or 
baryons (three quarks each with a different color). If two quarks separate 
far enough, the field energy increases and new quarks are created forming 
colorless hadrons. 

• Asymptotic freedom: At small distances the strength of the strong coupling 
is that low that quark and gluons behave as essentially free. This allows 
to use the perturbative approach in this regime, where a<j ^ 1. 



1.3 Deep inelastic scattering 

The scattering of electrons from protons, as illustrated in Figure II. 2[ has played 
a crucial role in the understanding of the proton structure. If the energy of the 
incoming electron {E) is low enough, the proton can be considered as a point 
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Figure 1.1: Summary of measurements of function of the energy scale Q, 

from [3j . 

charge (without structure). The differential cross section with respect to the 
solid angle of the scattered electron i^: 



where a (~ 1/137) is the fine structure constant, 9 is the angle at which the 
electron is scattered, E' is the outgoing electron energy and M the mass of the 
proton. E' is kinematically determined by 9. 

For higher energies of the incoming electrons, the interaction is sensitive to 
the proton structure, and the cross section becomes: 



Ki and K2 are functions that contain information on the proton structure and 
should be determined experimentally. Given that E' is kinematically determined 
by ^, Ki and K2 only dependent on one variable. 

^The mass of the electron is neglected in all formulas in this Section by assuming E » m. 




(1.3) 




(1.4) 
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k' = (E',k') 



k=(E,k) 




P' - (Po.P') 




Figure 1.3: Electron-proton deep inelastic scattering. 

Finally, for even higher electron energies, the proton breaks in a multi-hadron 
final state as illustrated in Figure 11.31 The cross section is then: 

Now p' is the sum of the momenta of the hadrons originating from the proton, 
and it is not constrained by p'^ = M^. Therefore, Wi and W2 are functions of 
two independent variables, E' and 9. Theoretically it is more convenient to use 
the Lorentz-invariant variables = —{k — k'Y and x = q'^/2qp, where p is the 
momenta of the incoming proton. 

The Parton Model describes the proton as built out of three point-like quarks 
('valence quarks') with spin 1/2, and interprets x as the fraction of the proton 
momentum carried by the quark. From the idea that at high the virtual 
photon interacts with a quark essentially free, Bjorken predicted that Wi and 
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W2 depend only on x at large (g^ > 1 GeV): 



MWi{q\x) Fi{x) 



(1.6) 




(1.7) 



According to the Parton Model: 



(1.8) 



where fi{x), called Parton Distribution Function (PDF), is the probability that 
the ith quark carries a fraction of the proton momentum x, and Qi is the electrical 
charge of the quark. Therefore, it is expected that 



but it was found experimentally that the result of this integral is 0.5. The rest of 
the proton momentum is carried by gluons. The introduction of gluons leads to 
a more complex description of the protons structure: quarks radiate gluons, and 
gluons produce qq pairs ('sea quarks') or radiate other gluons. Figure [Ol shows 
the PDFs of the valence quarks of the proton, the gluon, and the sea quarks. 
The valence quarks dominate at large x, whereas the gluon dominates at low x. 

The radiation of gluons results in a violation of the scaling behavior of Fi 
and F2, introducing a logarithmic dependence on g^, which is experimentally 
observed (see Figure [L6|) . The functional form of the PDFs can not be predicted 
from pQCD, but it is possible to predict their evolution with g^. 

The parton interactions at first order in are gluon radiation {q — )■ qg), 
gluon splitting {g — )• gg) and quark pair production {g — )■ qq). The probability 
that a parton of type p radiates a quark or gluon and becomes a parton of type 
p', carrying fraction y = x/z oi the momentum of parton p (see Figure [T3]) is 
given by the splitting functions: 




(1.9) 



+ 



y 



(1.10) 



y 



i-y 



3 y 



(1.11) 



PM = ^[y' + i^-yr] 



(1.12) 
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Figure 1.4: Example of PDFs of the valence quarks of the proton, the gluon, and 
the sea quarks as a function of x. 
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Figure 1.5: Diagrams at LO of the different parton interactions. 



= 3 



;i.i3) 



The evolution of the PDFs as a function of follow the DGLAP (Dokshitzer, 
Gribov, Lipatov, Altarelli and Parisi) equations |1]: 



(ilog(g2) 27r Ji 



X \ dz 
z J z 

X \ dz 



(1.14) 



;i.i5) 



d\og{q'^) 27r Ji \'—' ' ' 'z' ' ' ' ' ' z' j z 
The first equation describes the evolution of the quark PDF with due to 
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gluon radiation and quark pair production, whereas the second equation describes 
the change of the gluon PDF with due to gluon radiation and gluon sphtting. 
The equations assume massless partons and therefore are only valid for gluons 
and the light quarks (u, d and s). 



ZEUS 




1 10 lo' 10^ lo'' 7 iff' ^ 

e (GeV") 

Figure 1.6: Structure function F2 of the proton as measured by ZEUS, BCDMS, 
E665 and NMC experiments. 



1.4 Perturbative QCD 
1.4.1 The factorization theorem 

The QCD factorization theorem is a crucial concept of QCD, that states that 
cross sections in hadron-hadron interactions can be separated into a a hard par- 
tonic cross section (short-distance) component and a long-distance component. 
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described by universal PDFs: 

cr(Pi,P2) = J2j dxidx2fiixi, iJ,l)fj{x2, fijp) X aij{xi,X2,asilJ,l,fiR),q'^/lJ,%) 

(1.16) 

where Pi, P2 are the momenta of the interacting hadrons, the sum runs over all 
parton types, and (jjj is the partonic cross section of the incoming partons with 
hadron momenta fraction Xi, X2- Hr is the scale at which the renormalization is 
performed, and /^i? is an arbitrary parameter that separates the hard from the 
soft component. Both scales are typically chosen to be of the order of q^. 

Partonic cross sections in leading order (LO) calculations for jet production 
are O(a^), since they are based on 2 — )■ 2 parton interactions {gg gg,qg — >■ 
qg, qq — )■ qq), as shown in Figure [TT71 The dominant process is the gg scattering 
because of the larger color charge of the gluons. 

Next-to-leading-order (NLO) diagrams include contributions from gluon initial- 
or final-state radiation and loops on the diagrams already shown. The partonic 
cross sections at NLO reduce the dependence on the normalization and fac- 
torization scales, and are calculable using programs such as JETRAD [3] and 
NLOJET++ [6]. Predictions at higher orders are not yet available due to the 
large number of diagrams involved. 




Figure 1.7: Leading order diagrams for 2 —t- 2 parton interactions. 
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1.4.2 Parton Distribution Functions 

As already explained, perturbative QCD (pQCD) can predict the evolution of the 
PDFs [7J with respect to using the DGLAP equations, but not their functional 
form. Therefore, PDFs should be extracted from experimental data at a given 
= Ql- In particular, seven functions should be determined, one for the gluon 
and the others for each one of the light quarks and anti-quarks. Experimental 
data from a large variety of processes is used to constrain several aspects of the 
PDFs: measurements of Drell-Yan production, inclusive jet cross sections and 
W-asymmetry in pp collisions, and deep-inelastic e, n or u scattering. 

Typically, specific functional forms are postulated for the PDFs with a set 
of free parameters. These parameters are obtained optimizing the comparison 
between experimental data and predictions using the PDFs, for example by min- 
imizing a x^- The functional form assumed for several sets of PDFs is: 

f,{x,Ql)=x"^{l-xf'g,{x) (1.17) 

where and Pi are the free fit parameters and gi{x) is a function that tends 
to a constant in the limits a; — )■ and a; — ?■ 1. This functional form is mo- 
tivated by counting rules [8j and Regge theory |9j|, that suggest that fi{x) ~ 
(1 — x)^* when x — )■ 1 and fi{x) ~ x"* when x — )■ respectively. Both limits are 
approximate, and even if these theories predict the values of Pi and ai, they are 
taken as free fit parameters when computing the PDFs. This approach is used 
by three of the PDFs used in the analyses presented in this Thesis: CTEQ [TU] . 
MSTW HI] and HERA pjj PDFs. For example in the case of HERAPDFs, gi{x) 
is: 

gi{x) = 1 + eiX^/'^ + DiX + EiX^ (1.18) 

NNPDFs [13] follow a different approach, using neural networks as a parton 
parametrization. Neural networks are functional forms that can fit a large variety 
of functions. 

1.4.3 Uncertainties 

There are three main sources of uncertainties in the calculation of pQCD observ- 
ables: 

• The lack of knowledge of higher order terms neglected in the calculation. It 
is estimated by varying the renormalization scale, /ir, usually by a factor 



11 



Chapter 1. QCD at Hadron Colliders 



of two with respect to the default choice. The factorization scale, hf, 
is independently varied to evaluate the sensitivity to the choice of scale 
where the PDF evolution is separated from the partonic cross section. The 
envelope of the variation that these changes introduce in the observable is 
taken ClS db systematic uncertainty. 

• Uncertainties on parameters of the theory, like the as and the heavy quark 
masses, that are propagated into the observable. 

• PDFs have another uncertainty coming from the way the experimental data 
is used to determine the PDFs. This uncertainty is typically estimated 
using the Hessian method. If ao is the vector of the PDF parameters where 
X^(cio) is minimized, all parameters such that ~ Xo < T are considered 
acceptable fits, where T is a parameter called tolerance. PDF parameters 
are expressed in terms of an orthogonal basis, and variations along the 
positive and negative directions of each eigenvector ( ') such that 

— Xq = T are performed. The uncertainty in the observable F is: 



6T+ = /^max(r(a+) -r(ao),r(a-) -r(ao),0)2 (1.19) 



^F- = ^5]mm(F(a+) -F(ao),F(ar) -F(ao),0)2 (1.20) 

where F(a) is the observable computed using the PDFs with the parameters 
in vector a. NNPDF use a Monte Carlo approach to evaluate the uncer- 
tainties, in which the probability distribution in parameter space derives 
from a sample of MC replicas of the experimental data. Figure 11.81 shows 
the PDF of the gluon with its uncertainties obtained following different 
approaches. 



1.5 Monte Carlo simulation 

Complete pQCD calculations are always performed only up to a fixed order in 
as, but the enhanced soft-gluon radiation and coUinear configurations at higher 
orders can not be neglected. They are taken into account in the parton shower 
(PS) approximation, that sum the leading contributions of these topologies to all 
orders. Monte Carlo (MC) generator programs include the PS approximation, as 
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X 



Figure 1.8: PDF of the gluon as a function of x according to different PDF 
groups at = 2 GeV^. 

weU as models to reproduce non-perturbative effects, such as the hadronization 
of the partons to colorless hadrons and the underlying event (UE). 



1.5.1 Parton Shower 

The PS approximation describes successive parton emission from the partons in 
the hard interaction, as illustrated in Figure 11.91 The evolution of the showering 
is governed by DGLAP equations 11.141 and II. 15^ from which the Sudakov form 
factors [H] are derived for the numerical implementation of the parton shower. 
These factors represent the probability that a parton does not branch between 
an initial scale (tj) and a lower scale (t). Once a branching occurs at a scale 
ta, a — )■ 6c, subsequent branchings are derived from the scales and t^. They 
can be angle-, Q^- or p^-ordered. In the first case subsequent branchings have 
smaller opening angles than this between h and c, whereas in the second, parton 
emissions are produced in decreasing order of intrinsic pt- 

Successive branching stops at a cutoff scale, to, of the order of Aqcd, after 
producing a high-multiplicity partonic state. Since quark and gluons can not 
exist isolated, MC programs contain models for the hadronization of the partons 
into colorless hadrons. 
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Figure 1.9: Illustration of the parton shower from the outgoing partons of the 
hard interaction. 

1.5.2 Hadronization 

The hypothesis of local parton-hadron duality states that the momentum and 
quantum numbers of the hadrons follow those of the partons. This hypothesis 
is the general guide of all hadronization models, but do not give details on the 
formation of hadrons, that is described in models with parameters that are tuned 
to experimental data. There are two main models of hadron production. 

The string model [IB] describes the behavior of qq pairs using string dy- 
namics. The field between each qq pair is represented by a string with uniform 
energy per unit length. As the q and the q move apart from each other and thus 
the energy of the color field increases, the string connecting the two is tight- 
ened, until it breaks into a new q'q' pair. If the invariant mass of either of these 
string pieces is large enough, further breaks may occur. In the string model, the 
string break-up process is assumed to proceed until only on-mass-shell hadrons 
remain. In the simplest approach of baryon production, a diquark is treated 
just like an ordinary antiquark. A string can break either by quark-antiquark 
or antidiquark-diquark pair production, leading to three-quark states. There are 
more sophisticated models, but the formation of baryons is still poorly under- 
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stood. 

The cluster model [T7] is based on the confinement property of perturbative 
QCD, exploited to form color- neutral clusters. After the perturbative parton 
showering, all gluons are split into light quark-antiquark or diquark-antidiquark 
pairs. Color-singlet clusters are formed from the quarks and anti-quarks. The 
clusters thus formed are fragmented into two hadrons. If a cluster is too light to 
decay into two hadrons, it is taken to represent the lightest single hadron of its 
fiavor. Its mass is shifted to the appropriate value by an exchange of momenta 
with a neighboring cluster. If the cluster is too heavy, it decays into two clusters, 
that are further fragmented into hadrons. 

1.5.3 Underlying Event 

The UE comes from the partons that do not participate in the hard interaction. 
They contribute to the final state via their color- connect ion to the hard inter- 
action, and via extra parton-parton interactions. Its simulation is based on the 
eikonal model, that describes the underlying event activity as additional uncor- 
related partonic scatters. The number of interactions per event < n > depends 
on the impact parameter h. A small h value corresponds to a large overlap be- 
tween the two colliding hadrons, and therefore a higher probability for multiple 
interactions. For a given 6, the parton-parton cross section ahard is computed as 
a function of the transverse momentum in the center-of-mass frame of the scat- 
tering process pT- Since this cross section diverges as j)t — > a cut-off parameter 
p™" is introduced, where experimentally p™" ~ 2 GeV. < > is extracted from 
the ratio between the total hadron cross section and and the parton-parton cross 
section, < n >= (Jhard/cnd, and assumed to be Poisson-distributed. 

The UE models are tuned using experimental data, such as the jet shapes 
described in Chapters 4 and 5. 

1.5.4 Monte Carlo Generator Programs 
PYTHIA Monte Carlo 

The PYTHIA [18] MC event generator includes hard processes at LO, and uses 
the PS model for initial- and final-state radiation. The hadronization is per- 
formed using the string model. It includes an underlying event model to describe 
the interactions between the proton remnants. 
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The PYTHIA tunes DW [19] and Perugia2010 [20] use CTEQ5L PDFs, 
and both have been produced using Tevatron data. In the former the PS is 
Q^-ordered, whereas in the latter it is pr-ordered. 

In autumn 2009, the MRST LO* PDFs [21] were used in PYTHIA for the 
first time in ATLAS. This required to tune the PYTHIA model parameters, 
resulting in the MC09 [22j tune. It was done using Tevatron data, mainly from 
underlying event and minimum bias analyzes. The PS is p^-ordered. 

The PYTHIA-AMBTl [23] tune followed the MC09 one, and also uses 
MRST LO* PDFs and p^-ordered PS. It was derived using ATLAS data, in 
particular charged particle multiplicities in pp interactions at 0.9 and 7 TeV 
center-of-mass energy. 



HERWIG 

HERWIG [21] is a general-purpose MC event generator for hard processes in 
particle colliders. It uses an angular-ordered parton-shower for initial- and final- 
state QCD radiation, and a cluster model to reproduce the hadronization of the 
partons. The Fortran version of HERWIG is interfaced with JIMMY [25] to 
simulate multiple parton-parton interactions. 

HERWIG++ [26] is the C++ version of HERWIG, that is expected to replace 
the Fortran one at a given point. The underlying event is modeled inside the 
program, that therefore do not use JIMMY. 



ME + Parton Shower: Alpgen and Powheg 

Alpgen [27] is an event generator of multi-parton hard processes in hadronic col- 
lisions, that performs the calculation of the exact LO matrix elements for a large 
set of parton-level processes. It uses the ALPHA algorithm [28] to compute the 
matrix elements for large parton multiplicities in the final state. The advantage 
of this algorithm is that its complexity increases slower than the Feynman dia- 
grams approach when increasing the particles in the final state. Powheg [29] is 
a MC event generator that includes NLO matrix elements. Alpgen and Powheg 
contain an interface to both PYTHIA and HERWIG for the parton showering, 
the hadronization and the underlying event simulation. 
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1.6 Jet Algorithms 

Quarks and gluons from the hard scattering result on a collimated flow of par- 
ticles due to parton shower and hadronization. This collimated flow of particles 
is called jet. There are several jet deflnitions [30] with the main purpose of re- 
constructing jets with kinematics that reflect that of the initial parton. These 
deflnitions can be classifled in two main types of jet algorithms: cone algorithms 
and sequential recombination algorithms. 

1.6.1 Cone algorithms 

Typically, cone jet algorithms start by forming cones of radius R in the y—4> space 
around a list of seeds, that can be all particles in the flnal state or those above a 
given energy threshold. The center of the cone is recomputed from the particles 
inside by following either the snowmass or the four-momentum recombination. 
In the four-momenta recombination, the jet momenta is the sum of the momenta 
of its constituents: 

(E,p,,p„p,r*= j2iE,p.,py,p^y (1-21) 

const. 

whereas in the snowmass scheme, the jet is considered massless, its transverse 
energy is the sum of the transverse energy of its constituents and the jet [r], 0) 
are the average of the {rj, (p) of the constituents weighted by its transverse energy: 

= E't (1.22) 

const. 

(^,0r* = i T.iv,^yEfr (1.23) 

-t^T const. 

m?^^ = (1.24) 

A cone is formed from the new center and the process repeated until the 
particles inside the cone are no longer changed by further iterations. Usually 
the algorithm is allowed to form overlapping cones and then decides whether to 
merge or split them depending on the fraction of energy they share. 

This last step makes the cone algorithms collinear or infrared unsafe, and 
affects the deflnition of the parton-level jet cross section to all orders in pQCD. A 
jet algorithm is infrared safe if the addition of an extra particle with inflnitesimal 
energy do not change the jet conflguration in the final state. If the replacement 
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of a particle by two collinear particles (which momenta sum is equal to that of 
the original particle) do not change the jet configuration in the final state, the 
jet algorithm is collinear safe. 

In order to solve this, cone-based jet algorithms have been formulated such 
that they find all stable cones through some exact procedure, avoiding the use 
of seeds. These algorithms are very time-consuming from the computational 
point of view, which constitutes a disadvantage in high-multiplicity events such 
as those at the LHC 

1.6.2 Sequential recombination jet algorithms 

Sequential recombination jet algorithms cluster particles according to their rel- 
ative transverse momentum, instead of spacial separation. This is motivated by 
the parton shower evolution as described in Section 11.5.11 For all particles in the 
final state, the algorithm computes the following distances: 



where ku is the transverse momentum of particle i, Rij = y/Ay^~+~A^ between 
particles i and j, Ra. parameter of the algorithm that approximately controls the 
size of the jet, and p depends on the jet algorithm: p = 1 for the kt algorithm, 
p = for the Cambridge/ Aachen algorithm, andp = — 1 for the anti-fct algorithm. 
The distance diB is introduced in order to separate particles coming from the 
hard interaction than those coming from the interaction between remnants. The 
smallest distance is found, and if it is dij, particles i and j are combined into 
one single object. If instead it is diB, particle i is considered a jet an removed 
from the list. The distances are recalculated with the remaining objects, and 
the process repeated until no particle is left in the list. Jets are defined as those 
objects with pt above a given threshold. 

These algorithms are very convenient, mainly because they are infrared and 
collinear safe and computationally fast. In particular, the anti-kt algorithm [3T] 
produces jets with a conical structure in {y, 0), as illustrated in Figure [LT0| that 
facilitates dealing with pile-up. It is the default jet finding algorithm in the LHC 
experiments. 





l2p 
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Figure 1.10: A sample parton-level event clustered with the anti-/ct algorithm. 
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Chapter 2 



The ATLAS Detector at the 
Large Hadron Collider 



The analyses described in this Thesis are performed using proton-proton colhsion 
data produced by the Large Hadron CoUider (LHC) and collected by the ATLAS 
detector. In this Chapter, the LHC and the ATLAS detector are described, giving 
more emphasis to the elements that are relevant for the analyses. 



2.1 The Large Hadron Collider 

The LHC [32] is a superconducting accelerator built in a circular tunnel of 27 km 
in circumference that is located at CERN. The tunnel is situated between 45 to 
170 m underground, and straddles the Swiss and French borders on the outskirts 
of Geneva. Two counter rotating proton beams injected into the LHC from the 
SPS accelerator at 450 GeV are further accelerated to 3.5 TeV while moving 
around the LHC ring guided by magnets inside a continuous vacuum. During 
2010, the instantaneous luminosity was increased over time, with a maximum 
peak at 2 ■ 10^^ cm~^s~^, and the total integrated luminosity delivered by the 
LHC was of 48 pb~^ from which ATLAS recorded 45 pb~^ (see Figure [2?T]) . 

There are four main detectors placed along the accelerator line: ATLAS 
and CMS, that are general-purpose detectors, ALICE, dedicated to heavy-ions 
physics, and LHCb, dedicated to B-physics. 
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Day in 2010 Day in 2010 

Figure 2.1: Maximum instantaneous luminosity (left) and cumulative integrated lu- 
minosity (right) versus day delivered by the LHC and recorded by ATLAS for pp 
collisions at 7 TeV center-of-mass energy during stable beams in 2010. 

2.2 The ATLAS experiment 

The ATLAS detector [33] is an assembly of several sub-detectors arranged in 
consecutive layers around the beam shown in Figure 12.21 The main sub- 

detectors are the Inner Detector, the Calorimeters and the Muon System, that 
are described in the next Sections. ATLAS is 46 m long, 25 m wide and 25 m 
high, and weights 7000 t. 

The ATLAS coordinate system and its nomenclature will be used repeatedly 
throughout this Thesis, and is thus described here. The ATLAS reference system 
is a Cartesian right-handed coordinate system, with the nominal collision point 
at the origin. The anti-clockwise beam direction defines the positive 2;-axis, while 
the positive x-axis is defined as pointing from the collision point to the center of 
the LHC ring and the positive |/-axis pointing upwards. The azimuthal angle is 
measured around the beam axis, and the polar angle 9 is measured with respect 
to the 2;-axis. The pseudorapidity is defined as ?7 = — ln(tan(^/2)). The rapidity 
is defined as y = 0.5 x \n.[{E + pz) / {E ~Pz)], where E denotes the energy and pz 
is the component of the momentum along the beam direction. 

The ATLAS detector was designed to optimize the search for the Higgs boson 
and a large variety of physics phenomena at the TeV scale proposed by models 
beyond the Standard Model. The main requirements that follow from these goals 
are: 

• Given the high LHC luminosity, detectors require fast, radiation-hard elec- 
tronics and sensor elements. In addition, high detector granularity is 



22 



2.3. Inner Detector 



needed to handle the particle fluxes and to reduce the influence of overlap- 
ping events. 

• Large acceptance in pseudorapidity with almost full azimuthal angle cov- 
erage. 

• Good charged-particle momentum resolution and reconstruction efficiency 
in the inner tracker. 

• Very good electromagnetic calorimetry for electron and photon identifica- 
tion and measurements, complemented by full-coverage hadronic calorime- 
try for accurate jet and missing transverse energy measurements. 

• Good muon identification and momentum resolution over a wide range of 
momenta and the ability to determine unambiguously the charge of high 
Pt muons. 

• Highly efficient triggering on low transverse-momentum objects with suffi- 
cient background rejection is a prerequisite to achieve an acceptable trigger 
rate for most physics processes of interest. 



2.3 Inner Detector 

The Inner Detector (ID) was designed in order to perform high precision mea- 
surements with fine detector granularity in the very large track density events 
produced by the LHC The ID, that is ±3512 m long and 1150 mm in radius, 
is built out of three components, in increasing order of distance with respect to 
beam axis: the Pixel detector, the Semiconductor Tracker (SCT) and the Transi- 
tion Radiation Tracker (TRT). The precision tracking detectors (pixels and SCT) 
cover the region \ri\ < 2.5, and are segmented in r — and z, whereas the TRT 
cover the region |?7| < 2 and is only segmented in r — 0. The ID has around 87 
million readout channels, 80.4 millions in the pixel detector, 6.3 millions in the 
SCT and 351 thousand in the TRT. All three are immersed in a 2 T magnetic 
field generated by the central solenoid, which extends over a length of 5.3 m with 
a diameter of 2.5 m. 

The ID is used to reconstruct tracks and production and decay vertices, and 
provides a position resolution of 10, 17 and 130 jmv (Pixel, SCT, TRT) in the 
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Figure 2.2: View of the full ATLAS detector. 

r — (f) plane as well as 115 and 580 /ini (Pixel, SCT) in the r — z plane. The 
momentum resolution as a function of pt of the track is parametrized as: 

^ = Pl©P2XpT (2.1) 

Pt 

and the values Pi = 1.6 ±0.1% and P2 = (53 ±2) x 10~^ GeV^^ were determined 
using cosmic rays Extrapolation of the fit result yields to a momentum 

resolution of about 1.6% at low momenta and of about 50% at 1 TeV. 



2.4 Calorimeters 

The calorimeter systems of ATLAS, illustrated in Figure 1231 surround the Inner 
Detector system and cover the full 0-space and |?7| < 4.9, extending radially 
4.25 m. The calorimeter systems can be classified in electromagnetic calorime- 
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ters, designed for precision measurements of electrons and pliotons, and hadronic 
calorimeters, that collect the energy from hadrons. Calorimeter cells are pseudo- 
projective towards the interaction region in rj. The granularity of the electromag- 
netic calorimeter is typically 0.025 x 0.025 in lAr^l x |A<^|, whereas the hadronic 
calorimeters have granularity of 0.1 x 0.1 in most of the regions. The energy 
response of the calorimeter to single particles is discussed in the next Chapter. 



2.4.1 Liquid Argon Calorimeter 

The electromagnetic calorimeter is a lead-LAr detector with accordion-shaped 
kapton electrodes and lead absorber plates over its full coverage. The accor- 
dion geometry provides complete symmetry without azimuthal cracks. The 
calorimeter is divided into a barrel part {\r]\ < 1.475) and two end-cap compo- 
nents (1.375 < |?7| < 3.2), each housed in their own cryostat. Over the central 
region (|?7| < 2.5), the EM calorimeter is segmented in three layers in depth, 
whereas in the end-cap it is segmented in two sections in depth. Figure 12.41 
shows an sketch of a module of the LAr calorimeter. 

2.4.2 Hadronic calorimeters 

The Tile Calorimeter is placed directly outside the electromagnetic calorimeter 
envelope. Its barrel covers the region \ri\ < 1.0, and its two extended barrels the 



Tile barrel 



Tile extended barrel 




Figure 2.3: View of the calorimeter system. 
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Figure 2.4: Sketch of a module of the LAr calorimeter. 



range 0.8 < |?7| < 1.7. It is a sampling calorimeter using steel as the absorber 
and scintillating tiles as the active material, and it is segmented in depth in three 
layers, approximately 1.5, 4.1 and 1.8 interaction lengths (A) thick for the barrel 
and 1.5, 2.6, and 3.3 A for the extended barrel, as illustrated in Figure 1231 Two 
sides of the scintillating tiles are read out by wavelength shifting fibers into two 
separate photomultiplier tubes. 

The Hadronic End-cap Calorimeter (HEC) consists of two indepen- 
dent wheels per end-cap, located directly behind the end-cap electromagnetic 
calorimeter and sharing the same LAr cryostats. They cover the region 1.5 < 
|?7| < 3.2, and each wheel is divided into two segments in depth, for a total of 
four layers per end-cap. The HEC is a sampling calorimeter built out of copper 
plates intervealed with LAr gaps, that are the active medium. 

The Forward Calorimeter (FCal) is integrated into the end-cap cryostats, 
and it extends in |?7| from 3.1 to 4.9. It consists of three modules in each end-cap: 
the first, made of copper, is optimized for electromagnetic measurements, while 
the other two, made of tungsten, measure predominantly the energy of hadronic 
interactions. In all three modules LAr is the sensitive medium. 
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Figure 2.5: Barrel and extended barrel sections of the Tile Calorimeter. 
2.4.3 Calorimeter Topological Clusters 

Calorimeter clusters [35] are used as input to the jet finding algorithm in the stud- 
ies presented in this Thesis. They are built out of neighboring calorimeter cells 
with significant signal over noise. Therefore, the 3-D shape and the number of 
cells of clusters are not fixed. The noise is computed for each cell independently, 
and it is defined as the expected RMS of the electronic noise for the current gain 
and conditions plus the contribution from pileup added in quadrature. 

In order to make clusters, all cells with a signal to noise ratio above 4 are taken 
as seed cells. These cells are considered in descending order of signal to noise 
ratio, adding all neighboring cells to them forming the so called proto-clusters. 
Neighboring cells with signal to noise ratio between 2 and 4 are taken as seed 
cells in the next iteration. If a cell is adjacent to more than one proto-cluster 
and its signal to noise ratio is above 2 the proto-clusters are merged, whereas if 
it is smaller than 2 the cell is only added to the first proto-cluster. Once there 
are no more cells in the seed list, an splitting algorithm based on local maxima 
is run over the proto-clusters in order to separate clusters that are not isolated. 

Final clusters are treated as massless and their energy, at the electromag- 
netic (EM) scale, is the sum of the energies of the cells belonging to the cluster. 
The EM scale is the appropriate scale for the reconstruction of the energy de- 
posited by electrons or photons in the calorimeter. The EM energy scale was 
first determined using electron test-beam measurements in the barrel and endcap 



27 



Chapter 2. The ATLAS Detector at the Large Hadron Collider 



calorimeters [36] [37] . Muons from test-beams and produced in cosmic rays were 
used to validate the EM scale in the hadronic calorimeter [36] [38]. Recently, a 
recalibration of the electromagnetic calorimeters has been derived from Z — )■ ee 
events in pp collisions. The EM scale uncertainty is 1.5% in the EM LAr barrel 
calorimeter, and 3% in the Tile calorimeter. 



2.5 Muon System 

The calorimeter is surrounded by the Muon Spectrometer. The air-core toroid 
system, with a long barrel and two inserted end-cap magnets, generates a large 
magnetic field volume with strong bending power within a light and open struc- 
ture. Multiple-scattering effects are thereby minimized, and excellent muon 
momentum resolution is achieved with three stations of high-precision track- 
ing chambers. The muon instrumentation includes trigger chambers with very 
fast time response. 



2.6 Luminosity measurement 

The luminosity, £, of a pp collider that operates at a revolution frequency fr and 
rify bunches cross at the interaction point can be expressed as 

^ ^ nbfrnin2 2) 

where rii and n2 are the numbers of particles in the two colliding bunches and 
Hx and Tjy characterize the widths of the horizontal and vertical beam profiles, 
that are measured using van der Meer scans. The observed event rate is recorded 
while scanning the two beams across each other first in the horizontal (x), then 
in the vertical (y) direction. This measurement yields two bell-shaped curves, 
with the maximum rate at zero separation, from which one extracts the values 
of Sa; and S^. The luminosity at zero separation can then be computed using 
Equation 12.21 

ATLAS measures the luminosity [39] in inelastic interactions using different 
detectors and algorithms, all based on event-counting techniques, that is, on de- 
termining the fraction of bunch crossings (BCs) during which a specified detector 
registers an event satisfying a given selection requirement: 
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where /x is the average number of inelastic interactions per BC, atnei is the pp 
inelastic cross section, e is the efficiency for one inelastic pp collision to satisfy 
the event-selection criteria, and fiyis = e/i is the averaged number of visible 
(passing that event selection criteria) inelastic interactions per BC. The visible 
cross section = eamei is the calibration constant that relates the measurable 
quantity to the luminosity £ . Both e and a^s depend on the pseudorapidity 
distribution and particle composition of the collision products, and are therefore 
different for each luminosity detector and algorithm. 

In the limit /x^j^ <C 1, the average number of visible inelastic interactions per 
BC is given by the intuitive expression 

where N is the number of events passing the selection criteria that are observed 
during a given time interval, and Nbc is the number of bunch crossings in that 
same interval. When /i increases, the probability that two or more pp interactions 
occur in the same BC is not negligible, and is no longer linearly related to 
the raw event count A^. Instead must be calculated taking into account 
Poisson statistics, and in some cases, instrumental or pile-up related effects. 

a^s can be extracted from Equation 12.31 using the measured values of and 
£, computed with the van der Meer technique, that allows the determination of 
a^is without a priori knowledge of the inelastic pp cross section or of detector 
efficiencies. 



2.7 Trigger 

The trigger system has three different levels: Level-1 trigger (LVLl), Level-2 
trigger (LVL2) and the Event Filter (EF). Each trigger level refines the decisions 
made at the previous step and, when necessary, applies additional selection cri- 
teria. The rate of selected events is reduced to around 200 Hz for permanent 
storage trough the trigger chain, with an event size is approximately 1.3 Mbyte. 

The LI trigger searches for high transverse-momentum muons, electrons, pho- 
tons, jets, and r— leptons decaying into hadrons, as well as large missing and total 

^The LHC provides a collision rate of 40 MHz within full operational state. 
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transverse energy. Its selection is based on information from a subset of detectors, 
and reduces the event rate to about 75 kHz in less than 2.5 /is. In each event, the 
LI trigger defines one or more regions in 77 and from those regions within the 
detector where its selection process has identified interesting features. Data in 
these regions, called Regions-of-Intcrest (Rols), include information on the type 
of feature identified and the criteria passed. This information is subsequently 
used by the high-level trigger. 

The L2 selection is seeded by the Rol information provided by the LI trigger 
over a dedicated data path. L2 selections use, at full granularity and precision, all 
the available detector data within the Rols (approximately 2% of the total event 
data). The L2 menus are designed to reduce the trigger rate to approximately 
3.5 kHz, with an event processing time of about 40 ms, averaged over all events. 
The final stage of the event selection is carried out by the event filter, which 
reduces the event rate to roughly 200 Hz. Its selections are implemented using 
offline analysis procedures within an average event processing time of the order 
of 4 s. 

In the measurements presented in this Thesis, three different triggers have 
been used: the Minimum Bias Trigger Scintillators (MBTS), the central jet trig- 
ger, covering I77I < 3.2, and the forward jet trigger, spanning 3.1 < \r)\ < 4.9. 
The MBTS trigger requires at least one hit in the minimum bias scintillators 
located in front of the endcap cryostats, covering 2.09 < \r]\ < 3.84, and is the 
primary trigger used to select minimum-bias events in ATLAS. The central and 
forward jet triggers independently select data using several jet Et thresholds at 
the EM scale. More details on the trigger paths used to select events are given 
in next Chapter. 
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The jet cross section measurement is a fundamental test of QCD and it is sensi- 
tive to the presence of new physics. It also provides information on the parton 
distribution functions and the strong coupling. Since jet production has the 
largest cross section of all high pt processes at LHC, the related observables 
have the highest reach in energy at any given integrated luminosity. 

The first measurement of the inclusive jet cross section in ATLAS was done 
using an integrated luminosity of 17 nb~^ for jets with px > 60 GeV and \y\ < 
2.8, and was published in the EPJC [40j. In this Chapter, a measurement of 
the inclusive jet cross section for jets with pt > 20 GeV and \y\ < 4.4 [41j 
is presented, based on 37.3 ±1.2 pb^^ of integrated luminosity collected by 
ATLAS. The kinematic range of these two measurements is shown in Figure [3Jl 
as a function of the jet pt and \y\. In this Chapter, the MC samples used in 
the jet measurements are presented, followed by a detailed discussion of the 
jet reconstruction and calibration, the jet and event selection criteria, and the 
unfolding procedure to correct the measurement for detector effects. Finally, the 
measurement is compared to theoretical predictions. 

3.1 Monte Carlo simulation 

Samples of inclusive jet events in proton-proton collisions at -y/s = 7 TeV are 
used in the inclusive jet cross section analysis. PYTHIA 6.423 with the AMBTl 
tune MC samples are used to correct for detector effects, to compute the non- 
perturbative corrections applied to the NLO predictions, and to estimate part of 
the systematic uncertainties of the measurements. 
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Figure 3.1: Kinematic range of the inclusive jet cross section measured in this analy- 
sis [H] compared to that of the previous study reported in [ID] for jets identified using 
the anti-A;^ algorithm with R = 0.6. 

The MC generated samples are passed through a full simulation |12] of the 
ATLAS detector and trigger, based on GEANT4 that simulates the inter- 
actions of the generated particles with the detector material. The Quark Gluon 
String Precompound (QGSP) model 01] is used for the fragmentation of the 
nucleus, and the Bertini cascade (BERT) model [45j for the description of the 
interactions of the hadrons in the medium of the nucleus. Test-beam measure- 
ments for single pions have shown that these simulation settings best describe 
the response and resolution in the barrel [IB] and end-cap [H] calorimeters. 



3.2 Jet reconstruction and calibration 

In the analysis presented in this Chapter, jets are reconstructed using the anti-kt 
algorithm, with jet radius parameter R = 0.6 (results with R = 0.4 can be found 
in [II]). The jet finding algorithm is run over calorimeter clusters at the EM 
scale, and the four-momenta recombination is used. 

The EM scale needs to be further calibrated to account for calorimeter non- 
compensation (the energy response to hadrons is lower than the response to 
electrons of the same energy), dead material (inactive regions of the detector 
where energy is lost) and leakage (energy deposits from particles which shower 
is not fully contained in the calorimeter). Moreover, corrections are needed for 
particles inside the truth jet but not the calorimeter jet, and energy losses in 
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calorimeter clustering and jet reconstruction. In order to correct for all these 
effects, a jet energy scale (JES) correction as a function of the jet energy and 
pseudorapidity is applied to jets at the EM scale ^8] . 

The corrections are derived from Monte Carlo (MC) samples produced using 
PYTHIA 6.423 with the AMBTl tune. Truth jets are reconstructed from stable 
particles in the final state, excluding muons and neutrinos, that are taken into 
account when unfolding the measurement to the hadron level (see Section [37il) . 
Corrections are derived from pairs of calorimeter and truth jets matched within 
AR < 0.3, where AR = A/Arp~+~A0^. The energy response TZ is defined as 
the ratio of the calorimeter and the truth jet energy, TZ = ^ca«o^^trnt/t_ Both 
calorimeter and truth jets are required to be isolated, that is, only jets with no 
other jet with pt > 7 GeV at AR < 2.5R are accepted. The motivation for this 
requirement is that non-isolated jets have lower response than isolated jets, and 
the fraction of non-isolated jets in a given sample is physics dependent. The cut 
AR < 2.5R is chosen because anti-fct jets may extend a bit further than R. The 
average jet energy response, <7l >, and the average jet energy at the EM scale, 
< E'^"'^° >, are measured in bins of truth jet energy, E^"^^^^. A function is derived 
by fitting < 7^ > as a function of < E'^"'''" > with the following parametrization: 

n,aHb{E''^'°) = E aSniE'^^'ny (3-1) 

i 

where a, are free parameters of the fit. The JES correction is the inverse of this 
function, and it is applied to the jet four-momenta at the EM scale. 

Figure 13.21 shows the calorimeter EM response as a function of detector 
pseudo-rapidity, r/defl for several jet energies. The response has large fluctu- 
ations depending on the r]det region. Therefore, the procedure to derive the JES 
corrections is performed independently in different ridet bins. 

Figure 13.31 shows the average JES corrections as a function of pr in three rj 
regions. Corrections are larger in the central region (up to 2.1 at low px), and 
decrease as the px increases in all rj regions. 

An additional correction is required in cases where there are more than one 
proton-proton interaction in the same event (pile-up), because the jet energy 
includes other contributions apart from the one coming from the collision of 
interest. These corrections, called offset corrections [49j, have been derived by 
estimating the average extra-amount of energy per jet in data as a function of the 
number of primary vertices and rj. The offset corrections, shown in Figure 13. 4[ 

^rjdet is used for the jet pseudo-rapidity with respect to the center of the ATLAS detector. 
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Figure 3.2: Jet energy response at the EM scale as a function of the r]det, for several 
calibrated jet energy values. The intervals in ij^et used to evaluate the JES uncertainty 
are also shown. 



are subtracted from the jet energy at the EM scale before applying the JES 
correction. 

There are two corrections of the jet position that leave the jet energy unaf- 
fected. The first is needed because calorimeter clusters kinematics refer to the 
center of the ATLAS detector. In order to correct jets back to point to the pri- 
mary vertex, each jet constituent kinematics is recalculated using the vector from 
the primary vertex to the cluster centroid as its direction. The jet four-momenta 
is then re-set to the sum of the recalculated four-momenta of its constituents. 
This correction improves slightly (< 1%) the jet pr response. 

Finally, the jet i] is further corrected for the bias introduced by the variation 
of the jet energy response in the different parts of the calorimeter. This last cor- 
rection, added to the jet r], is smaller than 0.01 in most parts of the calorimeter, 
and goes up to 0.05 in the transition regions between the different calorimeters. 



3.3 Jet and event selection 

The analysis is based in the full 2010 data set of pp collisions at ^/s = 7 TeV 
collected by the ATLAS detector, corresponding to an integrated luminosity of 
37 pb~^. Events are considered when the trigger system, tracking detectors, 
calorimeters and magnets were operating at the nominal conditions. 

Events are then required to have at least one primary vertex with five or more 
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tative r] regions. Tlie correction is only shown over the valid kinematic range, E = 3.5 
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tracks pointing to it in order to remove beam-related backgrounds and cosmic 
rays. For events containing a jet with > 20 GeV, this requirement has an 
efficiency well above 99%. 

Jets with pt > 20 GeV and \y\ < 4.4 are selected, and required to pass 
quality criteria established to reject jets not coming from a proton-proton colli- 
sion. Fake jets may originate from calorimeter noise, mainly noise bursts in the 
hadronic endcap calorimeter electronics and coherent noise from the electromag- 
netic calorimeter, or from cosmic rays or beam background. For each source of 
fake jets a quality criteria has been designed by combining some of the following 
variables (for more details see [50]): 

• The fraction of the jet energy deposited in the electromagnetic calorimeter. 

• The maximum energy fraction in one calorimeter layer. 

• The energy fraction in the hadronic end-cap calorimeter. 

• Jet quality, defined as the fraction of LAr cells with a cell Q-factor greater 
than 4000. The cell Q-factor measures the difference between the measured 
pulse shape and the predicted pulse shape that is used to reconstruct the 
cell energy. 

• HEC jet quality, defined as the LArQuality except that is calculated only 
with the HEC calorimeter. 
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Figure 3.4: Jet offset at the EM scale shown as a function of pseudorapidity and the 
number of reconstructed primary vertices. 

• Negative energy in the jet, due to an offset correction higher than the jet 
energy. 

• Jet time computed as the energy squared cells mean time. 

• Minimum number of cells containing at least 90% of the jet energy. 

• Jet charged fraction, defined as the ratio of the sum of the pt of the tracks 
associated to a jet divided by the calibrated jet pt- 

A tag-and-probe method was use to determine the efficiency of the jet quality 
criteria, as shown in Figure 1331 for two representative rapidity regions. For jets 
with Pt > 100 GeV, the efficiency is above 99% in all rapidity regions. The 
maximum inefficiency, around 8%, is for central jets at low px- Jet cross sections 
are corrected for the inefficiency when it is larger than 1%. 

The MBTS trigger is used to select jets between 20 GeV and 60 GeV, whereas 
central and forward jet triggers are used to select jets with pt > 60 GeV. In the 
first runs (Period 1), the jet trigger was using only the LI decision to reject 
events, whereas starting from September (Period 2) the L2 decision was also 
used for rejection. The threshold in the L2 is typically 15 GeV larger than the 
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Figure 3.5: Efficiency for jet identification as a function of px for anti-fcj jets witli R 
= 0.6 for the rapidity regions |y| < 0.3 (left plot) and 2.8 < \y\ < 3.6 (right plot). The 
black circles indicate the efficiency measured in-situ using a tag-probe method, while 
the blue squares and the shaded band indicates the parametrized central value and the 
systematic uncertainty on the efficiency obtained by varying the tag jet selection. The 
turn-on is due to the fact that jet cleaning cuts are harder for low-p^^ jets. 



LI threshold in the same trigger chain, except for the LI 10 GeV threshold that 
corresponds to the L2 15 GeV threshold. As explained in Section I^Tl the trigger 
thresholds are at the EM scale. Figure 13.61 shows the trigger efficiency as a 
function of the trigger threshold and the calibrated jet pt in four representative 
\y\ regions. 

Given the increase of the instantaneous luminosity during 2010, prescales 
have been applied to jet triggers with low Et thresholds in order to control 
the trigger rate. For each jet pr and rapidity bin in this analysis, a trigger is 
used to select jets such that is fully efficient (> 99%) and that the prescale is 
as low as possible. The trigger efficiency for a given energy threshold, E^, is 
computed with a sample of the events passing a trigger with a lower threshold, 
E^ < E^. The efficiency is defined as the fraction of jets passing the trigger 
with E^ threshold with respect to the number of jets passing the trigger with 
E^ threshold. 

Tables 3.1 and 3.2 show the trigger used in each jet pt and \y\ bin. In the 
region \y\ < 2.8, the central jet trigger (L1_J, L2_J) is used. The trigger has 
lower efficiency in 1.2 < \y\ < 2.1 than in the rest of the central region. This is 
due to the fact that between rj of 1.3 and 1.5, there is a transition between the 
barrel and the end-cap of the liquid argon calorimeter. The electronic signals 
from these two parts of the calorimeter are not added up at the trigger level 
to integrate over the full jet energy deposition. The logical OR of central and 
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forward (L1_JF, L2_JF) jet triggers is used in the rapidity range 2.8 < \y\ < 3.6, 
since none of them are fully efficient alone. The forward jet trigger is used in the 
bin 3.6 < \y\ < AA. 
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Figure 3.6: Combined L1+L2 jet trigger efficiency as a function of reconstructed jet 
Pt for anti-let jets witli R = 0.6 in four representative |y| regions are sliown for different 
L2 trigger thresholds. The trigger thresholds are at the electromagnetic scale, while the 
jet pt is at the calibrated scale. The highest trigger chain does not apply a threshold 
at L2, so its LI threshold is listed. 
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Period 1 


Period 1 


Period 2 


Period 2 


Pt (GeV) 


y < 1.2 II y > 2.8 


1.2 < y < 2.1 


y <l.2\\ y > 2.8 


l.2< y < 2.1 


60-80 


L1_J5 


L1_J5 


L2_J15 


L2_J15 


80-110 


L1_J15 


L1_J5 


L2_J30 


L2_J15 


110-160 


L1_J30 


L1_J15 


L2_J45 


L2_J30 


160-210 


L1_J55 


L1_J30 


L2_J70 


L2_J45 


210-260 


L1_J75 


L1_J55 


L2_J90 


L2_J70 


260-310 


L1_J95 


L1_J75 


L1_J95 


L2_J90 


310-400 


L1_J115 


L1_J95 


L1_J115 


L1_J95 


> 400 


L1_J115 


L1_J115 


L1_J115 


L1_J115 



Table 3.1: Triggers used in the central region, \y\ < 2.8. 





Period 1 


Period 1 


Period 2 


Period 2 


PT (GeV) 


2.8 < 


y < 3.6 


3.6< y < 4.4 


2.8 < 


y < 3.6 


3.6 < y < 4.4 


60-80 


L1_J10 


L1_FJ10 


L1_J10 


not used 


L2_FJ25 


80-110 


L1_J10 


L1_FJ10 


L1_J30 


L2_J25 


1 L2_FJ25 


L2_FJ25 


110-1(30 


L1_J30 II L1_FJ30 


L1_.J55 


L2_.J15 II L2_F.J15 


L2_F.J 15 


160-210 


L1_J55 


L1_FJ55 


L1_J55 


L2_J45 


L2_FJ45 


L2_FJ70 


> 210 


L1_J55 


L1_FJ55 


L1_J55 


L2_J70 


L2_FJ70 


L2_FJ70 



Table 3.2: Triggers used in the forward region, 2.8 < |y| < 4.4. 
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3.4 Unfolding to the particle level 

The measured differential cross section with respect to the jet pr and y is defined 
as follows: 

da N^"^" 
dprdy CApxAy 

where N^'^^^ is the number of jets in the Ap^ and Ay bins, and £ is the integrated 
luminosity. 

This cross section is corrected for detector effects back to the particle level 
with a bin-by-bin unfolding procedure, using PYTHIA 6.423 with the AMBTl 
tune simulated samples. In order to improve the agreement between the MC 
and the data, the MC is reweighted in order to scale its prediction, that uses 
a modified LO PDF, to that of a NLO PDF. Corrections are derived in each 
Pt and y bin by computing the ratio between the inclusive jet pt distributions 
using truth jets (including muons and neutrinos) and using calorimeter simulated 
jets. These corrections are multiplied to equation 13. 2[ and have values between 
0.8 and 1.0 in most jet px regions, and around 0.6 at very low and very high jet 
Pt- 



3.5 Systematic Uncertainties 

In this section several sources of systematic uncertainty on the measurement are 
considered. More emphasis is given to the explanation of the JES uncertainty, 
since it is the main contribution to the total uncertainty on the inclusive jet cross 
section measurement. 

3.5.1 JES uncertainty 

The JES uncertainty comes mainly from the uncertainty on the single particle 
response, limitations in the detector knowledge (such as the amount of dead 
material), and the physics models and parameters (mainly fragmentation and 
underlying event) in the MC event generator used to derive the JES corrections. 
Seven rjdet bins have been selected to estimate the JES uncertainty, taking into 
account the position of the different calorimeters (see Figure 13. 2p . The JES 
uncertainty is first calculated in the region Irjdetl < 0.8 and then propagated 
forward in rapidity using a jet pr balance technique, explained later in this 
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section. The reason is that the detector is better known in the central region, and 
that test-beam measurements to estimate the uncertainty due to the calorimeter 
response to single particles were only performed in the range \r]det\ < 0.8. 

Uncertainty on the calorimeter response to single particles 

The calorimeter response to single hadrons has been studied comparing the 
calorimeter energy response to the momentum of isolated tracks |ST]. Results 
have been obtained from pp collisions at ^/s = 0.9 and 7 TeV, finding an agree- 
ment between data and MC within 2-5% for tracks with momentum between 0.5 
to 20 GeV. The calorimeter response to particles with momentum in the range 
between 20 GeV and 350 GeV has been derived from test-beam measurements, 
in which pion beams with momentum up to 350 GeV were injected into a slice 
of the ATLAS detector. For particles with larger momentum, an additional un- 
certainty of 10% on top of the 350 GeV measured uncertainty is applied in order 
to take into account longitudinal leakage and calorimeter non-linearities. The 
uncertainty of the single particle response in the calorimeter has been used to 
derive the related jet energy scale uncertainty. This is done using MC, and vary- 
ing randomly the energy of the particles of the jets separately, and according to 
the values of the single particle response uncertainty. The JES calibration con- 
stants are then computed with these modified jets and the results are compared 
to the nominal constants in order to derive the uncertainty. The uncertainty on 
the jet pt increases from 1.5% to 4% with increasing jet pt- 

Uncertainties due to the detector simulation 

Since the JES calibration is derived using MC samples, deviations in the descrip- 
tion of the detector in the MC introduce an uncertainty on the JES. This has 
been estimated by varying the material budget in the simulation. In the particle 
range of the E/p studies using isolated tracks only the dead material in the inner 
detector has been varied, since the measurement is already performed with the 
ATLAS detector. For other particles, also changes in the calorimeter material 
have been simulated. The jet response in these modified samples is compared to 
the nominal one in order to estimate the systematic uncertainty, that is always 
below 1%. 

As explained in Section I2.4.3[ clusters are build by adding neighboring cells 
with signal over noise above certain thresholds. Therefore, miss-modeling of 
the calorimeter cell noise in the MC would lead to different cluster formation 
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that eventually would affect the jet reconstruction. The impact of possible miss- 
modeling of the noise has been evaluated by varying the signal over noise thresh- 
olds in the cluster reconstruction. The difference of the jet energy response 
for jets reconstructed with the modified noise thresholds and the nominal jets 
has been taken as a systematic uncertainty. This goes up to 2% for jets with 
Pt < 45 GeV, and is negligible at higher pt- 

Uncertainties due to the event modelling in the Monte Carlo genera- 
tors 

The jet energy response obtained using PYTHIA with the AMBTl tune, used 
to derive the corrections, has been compared with those obtained using ALP- 
GEN+HERW1G+ JIMMY and PYTHIA with the Perugia2010 tune. This es- 
timates the uncertainty due to the different QCD description contained in the 
MCs, such as fragmentation or underlying event. It goes from 2% at low pt to 
1% for Pt > 100 GeV. 

Uncertainties due to the relative calibration in the forward regions 

As already mentioned, the JES uncertainties are evaluated for jets with \ri\ < 0.8. 
These uncertainties are transferred to the forward region by using a Pt balance 
technique in which the pr of forward jets is measured relative to that of central 
jets in dijet events. The difference of the pr balance of a central and a forward 
jet in data and in MC is taken as a systematic uncertainty, that is added in 
quadrature to the JES uncertainty in the central region. This is the largest 
source of Tinccrtainty at low px in the forward region, going up to around 12%. 
For jets with pt > 100 GeV it is always smaller than 2%. 

Uncertainty of the method to derive the JES calibration constants 

There is an uncertainty associated to the method used to derive the JES calibra- 
tion constants. It comes mainly from the assumption that every jet needs the 
same average compensation and from the fact that the same correction factor is 
used for the energy and the pt, whereas the jet mass is not always zero. Two sta- 
tistical independent MC samples were used to derive the uncertainty associated 
to the method itself. One was used to derive the corrections, that were applied 
to the calorimeter jets of the other. Deviations of the jet response from unity 
in this second sample are taken as systematic uncertainties, that are generally 



43 



Chapter 3. Inclusive Jet Cross Section 



below 1%, apart from the central region, where the uncertainty amounts up to 
2% at low pt- 

Total uncertainty on the jet pt 

In Figure the fractional JES uncertainty and its components are shown as a 
function of Pt in three representative rj regions. The fractional JES uncertainty 
is below 4% for px > 80 GeV in all t] regions, it is smaller than 5% for all pt in 
the central region, and amounts up to 13% at low pt in the most forward region. 

Uncertainty due to multiple interactions 

A final contribution to account for the uncertainty due to pile-up is added in 
quadrature in events with more than one primary vertex. This uncertainty is 
estimated mainly from the variation of the offset corrections derived using tracks 
instead of calorimeter towers. It is shown in Figure 13.81 for events with two 
primary vertices. In this case, the uncertainty due to pile-up for low pt jets is 
around 1% in the central region and almost 2.5% in the most forward region. 
The uncertainty is always below 1% for jets with px > 60 GeV. In the case of 
three primary vertices, the pile-up uncertainty is approximately twice that of 
two primary vertices, and with four primary vertices the uncertainty for central, 
endcap and forward jets is less than 3%, 6% and 8%, respectively. The relative 
uncertainty due to pile-up for events with 5 additional interactions becomes less 
than 1% for all jets with pt > 200 GeV. 

Uncertainty on the measurement due to the JES uncertainty 

The JES uncertainty is the largest contribution to the uncertainty of the inclusive 
jet cross section measurement, that goes from around +30%/-20% in the central 
region to approximately +80%/-50% in the most forward region. 

3.5.2 Other sources of systematic uncertainties 

• The unfolding factors have been recomputed varying the jet pt in the MC 

in order to take into account the jet energy and angular resolution, and the 
differences on the cross section shape in data and in MC. These unfolding 
factors have been compared to estimate the systematic uncertainty, that is 
typically between 2% and 5%, except for the lowest pt bin where it goes 
up to 20%. 
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• Matching of track jets to calorimeter jets has been performed in both data 
and Monte Carlo to estimate the modehng of the calorimeter jet recon- 
struction efficiency in the MC. The difference is smaller than 2% (3%) for 
jets with pt > 20 GeV (30 GeV), and is taken as a systematic uncertainty. 

• In this analysis, jets are selected with triggers that are 99% efficient. A 
conservative 1% uncertainty overall has been considered to account for the 
trigger inefficiency. 

• The systematic uncertainty on the efficiency of the jet quality criteria 

is taken as a systematic uncertainty on the cross section, and it is always 
below 1%. 

• Finally, the uncertainty on the luminosity measurement is 3%. 
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Figure 3.7: Fractional JES uncertainty as a function of px for jets in different rj 
regions. Tlie JES uncertainty for regions \7]\ > 0.8 contains tlie contribution from the 
rj inter-calibration between central and forward jets in data and Monte Carlo added 
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3.6 Results 

The inclusive jet cross section unfolded to the hadron level is shown in Figure 13^ 
for different \y\ regions up to \y\ < 4.4. It is measured from 20 GeV to 1.5 TeV 
and spans up to ten orders of magnitude. The data are compared to NLO predic- 
tions corrected for non-perturbative effects. The NLO predictions are obtained 
from NLOJET-I-+4.1.2 program with CTEQ 6.6 PDFs, and the corrections for 
non-perturbative effects are derived using samples produced with PYTHIA 6.423 
with the AMBTl tune and applied to NLO predictions. These corrections are 
obtained bin-by-bin, comparing the cross section with and without hadroniza- 
tion and underlying event. The correction is dominated by the underlying event 
at low pt (1-5 at 20 GeV), and tends to 1 as the px increases. The uncer- 
tainty on the NLO predictions is obtained adding in quadrature uncertainties 
from the PDFs, the choice of factorization and renormalization scales, and the 
value of the strong coupling constant. It is combined with the uncertainty on the 
non-perturbative corrections, estimated from the maximum difference of the cor- 
rections with AMBTl with respect to other PYTHIA tunes and to HERWIG++. 
The total uncertainty in the predictions is typically around 20% at low and high 
Pt and 10% at intermediate pt values. 

Figure 13.101 shows the ratio of the measured cross section in data and the 
theory prediction, that are in agreement within uncertainties. The cross section 
in data is lower than the cross section predicted by NLO predictions in the 
forward region and at high px- 

In Figure [3?TT| data are compared to results using CTEQ 6.6, MTSW 2008, 
NNPDF 2.1, and HERAPDF 1.5. All cross sections are normalized to that 
obtained with CTEQ 6.6. Predictions using MTSW 2008, NNPDF 2.1, and 
HERAPDF 1.5 are closer to data than those using CTEQ 6.6, but all have the 
tendency to produce higher cross sections than data at high pt in the forward 
region. This may be due to the fact that PDFs are currently poorly constrained 
in this kinematic region. 

Finally, Figure 13. 121 shows the comparison of data and POWHEG predictions 
using MTSW 2008, both normalized to the NLO with MTSW 2008 cross section 
shown in Figure ETI POWHEG is interfaced with either PYTHIA or HERWIG 
for the parton showering, the hadronization and the underlying event modeling. 
POWHEG predictions agree with data within uncertainties, but tend to produce 
larger cross section at low px and, mainly POWHEG+HERWIG, smaller cross 
sections at high pt- POWHEG gives different results depending on whether 
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PYTHIA or HERWIG is used for the showering, which may be taken as an indi- 
cation of the level of uncertainty related to the parton shower implementation. 
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Figure 3.9: Inclusive jet cross section as a function of pT in different regions of |y| 
for jets identied using the anti-A;j algorithm with R = 0.6. For convenience, the cross 
sections are multiplied by the factors indicated in the legend. The data are compared to 
NLO QCD calculations to which non-perturbative corrections have been applied. The 
error bars indicate the statistical uncertainty on the measurement, and the shaded band 
indicates the quadratic sum of the systematic uncertainties. There is an additional 
overall uncertainty of 3% due to the luminosity measurement that is not shown. 
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Figure 3.11: Inclusive jet cross section as a function of pT in different regions of |y| 
for jets identied using the anti-A:f algorithm with R = 0.6. The theoretical error bands 
obtained by using different PDF sets (CTEQ 6.6, MSTW 2008, NNPDF 2.1, HERA 
1.5) are shown. The data points and the error bands are normalized to the theoretical 
estimates obtained by using the CTEQ 6.6 PDF set. 
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Figure 3.12: Inclusive jet cross section as a function of pT in different regions of \y\ 
for jets identied using the anti-Zc^ algorithm with R = 0.6. The ratios of the POWHEG 
predictions showered by PYTHIA and HERWIG and the data to the NLO predictions 
corrected for the non-perturbative effects is shown. The ratio shows only the statistical 
uncertainty on the POWHEG prediction. The total systematic uncertainties on the 
theory and measurement are indicated. 
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One of the fundamental elements of jet measurements is the proper understanding 
of the energy flow around the jet core and the validation of the QCD description 
contained in the event generators, such as parton shower cascades, and the frag- 
mentation and underlying event models. Jet shapes [52] observables are sensitive 
to these phenomena and thus very adequate to this purpose. They have been 
measured in pp [S3], e^p [51], and e~^e~ [SB] colliders. 

Before the first pp collisions delivered by the LHC, the sensitivity of jet shapes 
to the UE and pile-up in ATLAS was shown using MC simulation at 10 TeV [56] 
(see Appendix . Jet shapes studies at detector level were performed with the 
very first ATLAS data at a center of mass energy of 900 GeV [SZj, [5B] (see 
Appendix [B]) . 

Results at calorimeter level using 7 TeV data were used in the first measure- 
ment of the inclusive jet cross section in ATLAS [ID] based on 17 nb~^ to validate 
the QCD description in the event generators and to consolidate the confidence 
in the corrections and calibrations applied. Jet shapes at calorimeter level in the 
forward region {\y\ > 2.8) have been computed with the same purpose for the 
inclusive jet cross section measurement using the full 2010 dataset [H], corre- 
sponding to 37 pb^^. 

In this Chapter, measurements of differential and integrated jet shapes in 
proton-proton collisions at y/s = 7 TeV are presented [SH], using data collected 
by the ATLAS experiment corresponding to 3 pb~^ of total integrated luminos- 
ity. The definition of jet shapes observables and the event selection criteria are 
explained, followed by a description of Monte Carlo samples used and a detailed 
discussion on the corrections for detector effects and systematic uncertainties. Fi- 
nally, results are presented and compared to predictions from several MC event 
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generators. 



4.1 Jet shape definition 

The differential jet shape p(r) as a function of the distance r = ^/Ay^~+~A(p 
to the jet axis is defined as the average fraction of the jet pt that hes inside an 
annulus of inner radius r — Ar/2 and outer radius r + Ar/2 around the jet axis: 



11 ^prir - Ar/2,r + Ar/2) ^ , 
or = y ^ ' , 2 <r < R- Ar 2, 



(4.1) 



where PTij'i,r2) denotes the scalar sum of the px of the jet constituents in the 
annulus between radius ri and r2, A^-''^* is the number of jets, and R = 0.6 and 
Ar = 0.1 are used. The points from the differential jet shape at different r values 
are correlated since, by definition, Y^o pi^^) = 1- Alternatively, the integrated 
jet shape \E'(r) is defined as the average fraction of the jet px that lies inside a 
cone of radius r concentric with the jet cone: 

^(^) = ^yM^ 0<r<i?, (4.2) 

where, by definition, \E'(r = R) = 1, and the points at different r values are 
strongly correlated. This geometrical definition of the shape of the jet, based 
on concentric cones around the jet axis, is particularly adequate for jets recon- 
structed using the anti-fcf algorithm. 

In this analysis, the jet finding algorithm is run over uncorrected energy clus- 
ters or towers (calorimeter level). Inner Detector tracks, or final-state particles in 
the MC generated events (particle level). The jet shape measurements are per- 
formed in different regions of jet pt and \y\, and a minimum of 100 jets in data 
are required in each region to limit the statistical fluctuations on the measured 
values. 



4.2 Event selection 

The event and jet selection criteria, close to that of the inclusive jet cross section 
measurement detailed in Section [331 ^-I'e: 



54 



4.3. Monte Carlo simulation 



• Detector components relevant for this analysis, such as the calorimeter and 
the Inner Detector, operating at the nominal conditions. 

• Only data from runs until the end of August are used, in order to avoid 
bunch train pile-up 0. This corresponds to 3 pb~^ of total integrated lumi- 
nosity collected by the ATLAS experiment. 

• The trigger is required to be fully efficient in each pt and y bin. 

• The events are required to have one and only one reconstructed primary 
vertex with a z-position within 10 cm of the origin of the coordinate sys- 
tem, which suppresses pile-up contributions from multiple proton-proton 
interactions in the same bunch crossing, beam-related backgrounds and 
cosmic rays. The tight cut on the z-position is performed in order to keep 
the projective geometry of the calorimeter. 

• Only jets with corrected transverse momentum pt > 30 GeV and rapidity 
\y\ < 2.8 are considered. 

4.3 Monte Carlo simulation 

Monte Carlo simulated samples are used to correct the jet shapes for detector 
effects back to the hadron level, and to estimate part of the systematic uncertain- 
ties. PYTHIA 6.4.21, HERWIG++ 2.4.2 and ALPGEN 2.13 event generators are 
used to produce inclusive jet events in proton-proton collisions at a/s = 7 TeV. 
In the case of PYTHIA, samples are generated using three tuned sets of param- 
eters denoted as ATLAS-MC09, DW, and Perugia2010. In addition, a special 
PYTHIA-Perugia2010 sample without UE contributions is generated using the 
RIVET package [60]. HERWIG++ and PYTHIA-MC09 samples are generated 
with MRST2007LO* [II] [21] PDFs, PYTHIA-Perugia2010 and PYTHIA-DW 
with CTEQ5L [6T] PDFs, and ALPGEN with CTEQ61L [62j PDFs. 

4.4 Jet shapes at calorimeter level 

The measured differential (integrated) jet shapes, as determined by using calorime- 
ter clusters, are shown in Figures 14.11 and 14.21 (Figures 14.31 and 14. 4p . The data 

^ Bunch train pile- up refers to the energy deposit coming from pp colhsions in a different 
bunch crossing than the one of the event under consideration. 
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present a prominent peak at low r, which indicates that the majority of the jet 
momentum is concentrated around the jet axis. As expected, jets get narrower 
as the pt increases. The data are compared to several MC predictions, being 
PYTHIA-Perugia2010 the one that describes the data best. A reasonable agree- 
ment between jet shapes data and PYTHIA-Perugia2010 has been also observed 
in the different \y\ regions considered in this analysis. Therefore, PYTHIA- 
Perugia2010 is used to correct both differential and integrated jet shapes for 
detector effects back to the particle level. 

4.5 Correction for detector effects 

The correction for detector effects is done using a bin-by-bin procedure that also 
accounts for the efficiency of the selection criteria and of the jet reconstruction 
in the calorimeter. Here, the method is described in detail for the differential 
case. A similar procedure is employed to correct independently the integrated 
measurements. The correction factors U{r,pT,\y\) are computed separately in 
each jet pt and \y\ region. They are defined as the ratio between the jet shapes 
at the particle level p{r)P^^, obtained using particle-level jets in the kinematic 
range under consideration, and the reconstructed jet shapes at the calorimeter 
level after the selection criteria are applied and using calorimeter-level 

jets in the given pt and \y\ range. 

The correction factors U{r,pT, \y\) = p{.f)mc I Pi.''^)mc present a moderate pt 
and \y\ dependence and vary between 0.95 and 1.1 as r increases. For the inte- 
grated jet shapes, the correction factors differ from unity by less than 5%. This 
is shown in Figures 14.51 and 14.61 (14.71 and 14.81) for the differential (integrated) 
jet shapes. The corrected jet shape measurements in each pf and \y\ region are 
computed by multiplying bin-by-bin the measured uncorrected jet shapes in data 
by the corresponding correction factors. 
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Figure 4.1: The measured uncorrected differential jet sliape using calorimeter clusters 
for jets ■with \y\ < 2.8 and 30 GeV < px < 210 GeV is shown in different pT regions. 
The predictions of PYTHIA-Perugia2010 (sohd lines), HERWIG++ (dashed hues), 
ALPGEN interfaced with HERWIG and JIMMY (dotted hues), and PYTHIA-MC09 
(dashed-dotted lines) are shown for comparison. Only statistical uncertainties are 
considered. 



57 



Chapter 4. Jet Shapes 



10-^ 



Data 

- Pythia Perugia2010 

■ He™ig++ 

■ Alpgen 
Pythia MC09 



anti-kjets R = 0.6 

21 QeV < p^" < 260 QeV 

I I < 2.8 

Calorimeter 

(oniy statistical uncertainties) 



0.1 



0.2 



0.3 



0.4 



0.5 



0.6 



10-^ 



anti-I^Jets R = 0.6 
260QeV<p^"<310QeV 
l/'l<2-8 
Caiorimeter 

(oniy statistical uncertainties) 



Data 

- Pyttiia Perugia2010 

■ Herwig++ 

■ Alpgen 
Pyttiia MC09 



0.1 



0.2 



0.3 



0.4 



0.5 



0.6 



10' 



' Data 

— Pyttiia Perugia2010 

- - Herwig++ 
Alpgen 

■ - Pyttiia MG09 



anti-k, jets R = 0.6 

31 GeV < p^" < 400 QeV 

I /' I < 2.8 

Calorimeter 

(only statistical uncertainties) 



0.1 



0.2 



0.3 



0.4 



0.5 



0.6 



10-^ 



' Data 

— Pyttiia Perugia2010 

- - Herwig++ 
Alpgen 

■ - Pyttiia MC09 



anti-k Jets R = 0.6 
400 QeV < p^" < 500 QeV 
I I < 2.8 
Calorimeter 

(only statistical uncertainties) 



0.1 



0.2 



0.3 



0.4 



0.5 



0.6 



1 - 



10-' 



Data 

■ Pyttiia Perugia2010 

■ Henfl/ig++ 
Alpgen 
Pyttiia MC09 



anti-k, jets R = 0.6 
500 QeV < < 600 QeV 
I I < 2.8 
Calorimeter 

(only statistical uncertainties) 



4^ 



0.1 



0.2 



0.3 



0.4 



0.5 



0.6 



Figure 4.2: The measured uncorrected differential jet shape using calorimeter clusters 
for jets with \y\ < 2.8 and 210 GeV < px < 600 GeV is shown in different pT regions. 
The predictions of PYTHIA-Perugia2010 (soUd lines), HERWIG++ (dashed hues), 
ALPGEN interfaced with HERWIG and JIMMY (dotted hues), and PYTHIA-MC09 
(dashed-dotted lines) are shown for comparison. Only statistical uncertainties are 
considered. 
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Figure 4.3: The measured uncorrected integrated jet shape using calorimeter clusters 
for jets with \y\ < 2.8 and 30 GeV < px < 210 GeV is shown in different pT regions. 
The predictions of PYTHIA-Perugia2010 (sohd lines), HERWIG++ (dashed hues), 
ALPGEN interfaced with HERWIG and JIMMY (dotted hues), and PYTHIA-MC09 
(dashed-dotted lines) are shown for comparison. Only statistical uncertainties are 
considered. 
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Figure 4.4: The measured uncorrected integrated jet shape using calorimeter clusters 
for jets with \y\ < 2.8 and 210 GeV < pT < 600 GeV is shown in different pT regions. 
The predictions of PYTHIA-Perugia2010 (sohd lines), HERWIG++ (dashed hues), 
ALPGEN interfaced with HERWIG and JIMMY (dotted hues), and PYTHIA-MC09 
(dashed-dotted lines) are shown for comparison. Only statistical uncertainties are 
considered. 
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Figure 4.5: Correction factors applied to the measured differential jet shapes to correct 
the measurements for detector effects for jets with \y\ < 2.8 and 30 GeV < pT < 
210 GeV. 
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Figure 4.6: Correction factors applied to tire measured differential jet shapes to correct 
the measurements for detector effects for jets with \y\ < 2.8 and 210 GeV < pT < 
600 GeV 
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Figure 4.7: Correction factors applied to tlie measured integrated jet shapes to correct 
the measurements for detector effects for jets with \y\ < 2.8 and 30 GeV < pT < 
210 GeV. 
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Figure 4.8: Correction factors applied to tlie measured integrated jet shapes to correct 
the measurements for detector effects for jets with \y\ < 2.8 and 210 GeV < pT < 
600 GeV 
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4.6 Systematic uncertainties 

A detailed study of systematic uncertainties on the measured differential and 
integrated jet shapes has been performed. The impact on the differential mea- 
surements is described here. 

• The absolute energy scale of the individual clusters belonging to the jet 
is varied in the data according to studies using isolated tracks |40j, which 
parametrize the uncertainty on the calorimeter cluster energy as a function 
of pt and rj of the cluster. This variation is performed by multiplying the 
cluster energies by 1 ± a x (1 + for up (down) variation, with 
a = 5% if \v\ < 3.2 (a = 10% if |r/| > 3.2), b = 1.5, and p^"***^^ in GeV. 
The jet position is then recalculated according to the new cluster four- 
momenta. The maximum variation on the differential jet shapes is taken 
as a systematic uncertainty, that varies typically between 3% to 15% as 
r increases, and constitutes the dominant systematic uncertainty in this 
analysis (see Figures and H. lOp . 

• The systematic uncertainty on the measured jet shapes arising from the de- 
tails of the model used to simulate calorimeter showers in the MC events is 
studied. A different simulated sample is considered, where the FRITIOF |63] 
plus BERT showering model is employed instead of the QGSP plus BERT 
model. FRITOF-I-BERT provides the second best description of the test- 
beam results after QGSP+BERT. This introduces an uncertainty on 
the measured differential jet shapes that varies between 1% to 4%, and is 
approximately independent of pt and \y\ (see Figures 14.111 and l4.12p . 

• The measured jet pt is varied by 2% to 8%, depending on px and \y\, to 
account for the remaining uncertainty on the absolute jet energy scale |i4D] . 
after removing contributions already accounted for and related to the en- 
ergy of the single clusters and the calorimeter shower modeling, as discussed 
above. This introduces an uncertainty of about 3% to 5% in the measured 
differential jet shapes (see Figures and 14. lOp . 

• The 14% uncertainty on the jet energy resolution [40] translates into a 
smaller than 2% effect on the measured differential jet shapes (see Fig- 
ures SIS] and SH]) . 

• The correction factors are recomputed using HERWIG++, which imple- 
ments different parton shower, fragmentation and UE models than PYTHIA, 
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and compared to PYTHIA-Perugia2010 (see Figures 14.151 and 14.161) . In 
addition, the correction factors are also computed using ALPGEN and 
PYTHIA-DW for px < HO GeV, where these MC samples provide a rea- 
sonable description of the uncorrected shapes in the data. The results from 
HERWIG++ encompass the variations obtained using all the above gen- 
erators and are conservatively adopted in all pr and \y\ ranges to compute 
systematic uncertainties on the differential jet shapes. These uncertainties 
increase between 2% and 10% with increasing r. 

• Two statistically independent samples are used to test the closure of the 
correction for detector effects procedure. One was used to compute the 
shapes using calorimeter clusters, while the other was used to correct these 
shapes for detector effects. In Figures 14.171 and 14.181 corrected results from 
the first sample are compared to results at particle level from the second 
one. The ratio between the two deviates from unity about 1%, value that 
is included as uncertainty on the differential measurements. 

• No significant dependence on instantaneous luminosity is observed in the 
measured jet shapes, indicating that residual pile-up contributions are neg- 
ligible after selecting events with only one reconstructed primary vertex. 

• It was verified that the presence of dead calorimeter regions in the data 
does not affect the measured jet shapes, since only affected a small fraction 
of jets. 

The different systematic uncertainties are added in quadrature to the statisti- 
cal uncertainty to obtain the final result. The total systematic uncertainty for 
differential jet shapes decreases with increasing pt and varies typically between 
3% and 10% (10% and 20%) at r = 0.05 (r = 0.55) (see Figures gH and |i;20D . 
The total uncertainty is dominated by the systematic component, except at very 
large pt where the measurements are still statistically limited. In the case of the 
integrated measurements, the total systematic uncertainty varies between 10% 
and 2% (4% and 1%) at r = 0.1 (r = 0.3) as pr increases, and vanishes as r 
approaches the edge of the jet cone (see Figures 14.211 and l 4.22p . 
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Figure 4.9: Systematic uncertainty on the differential jet shape related to the absolute 
energy scale uncertainties on clusters and jets, for jets with |y| < 2.8 and 30 GeV < 
PT < 210 GeV. 
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Figure 4.10: Systematic uncertainty on tlie differential jet shape related to the ab- 
solute energy scale uncertainties on clusters and jets, for jets with |y| < 2.8 and 
210 GeV <pT< 600 GeV. 
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Figure 4.11: Systematic uncertainty on the differential jet shape related to the 
calorimeter showering model, for jets with \y\ < 2.8 and 30 GeV < px < 210 GeV. 
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Figure 4.12: Systematic uncertainty on the differential jet shape related to the 
calorimeter showering model, for jets with \y\ < 2.8 and 210 GeV < px < 600 GeV. 
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Figure 4.13: Systematic uncertainty on the differential jet shape related to the jet pT 
resolution, for jets with \y\ < 2.8 and 30 GeV < pr < 210 GeV. 
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Figure 4.14: Systematic uncertainty on the differential jet shape related to the jet pT 
resolution, for jets with \y\ < 2.8 and 210 GeV < pr < 600 GeV. 
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Figure 4.15: Systematic uncertainty on the differential jet shape related to the cor- 
rection for detector effects with different physics models assumptions, for jets with 
\y\ < 2.8 and 30 GeV < pr < 210 GeV. 
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Figure 4.16: Systematic uncertainty on the differential jet shape related to the cor- 
rection for detector effects with different physics models assumptions, for jets with 
\y\ < 2.8 and 210 GeV < pr < 600 GeV. 
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Figure 4.17: Systematic uncertainty on the differential jet shape related to the non- 
closure of the correction for detector effects procedure, for jets with \y\ < 2.8 and 
30 GeV <PT< 210 GeV. 
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Figure 4.18: Systematic uncertainty on the differential jet shape related to the non- 
closure of the correction for detector effects procedure, for jets with \y\ < 2.8 and 
210 GeV <PT < 600 GeV. 
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Figure 4.19: Summary of systematic uncertainties for the differential jet shape mea- 
surements for jets with \y\ < 2.8 and 30 GeV < px < 210 GeV. 
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Figure 4.20: Summary of systematic uncertainties for the differential jet shape mea- 
surements for jets with \y\ < 2.8 and 210 GeV < pT < 600 GeV. 
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Figure 4.21: Summary of systematic uncertainties for the integrated jet shape mea- 
surements for jets with \y\ < 2.8 and 30 GeV < px < 210 GeV. 
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Figure 4.22: Summary of systematic uncertainties for the integrated jet shape mea- 
surements for jets with \y\ < 2.8 and 210 GeV < pr < 600 GeV. 
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4.6.1 Cross-checks using other detector objects 

The jet shape analysis is also performed using tracks inside the cone of jets 
reconstructed using calorimeter clusters. The inner detector tracks are required 
to pass selection criteria derived in the analysis of charged particle multiplicities 
in ATLAS [M]: 

• Pt> 100 MeV and \r]\ < 2.5 

• a hit in the first layer (layer-0) of the Pixel detector if one is expected 
(meaning that the layer-0 module is not out of the configuration and the 
track is not close to the module edge) 

• a minimum of one Pixel hit in any of the 3 layers of this detector 

• at least two (if pr > 100 MeV), four (if pr > 200 MeV) or six (if pr > 
300 MeV) SCT hits 

• transverse and longitudinal impact parameters calculated with respect to 
the event primary vertex |(iOpv| < 1.5 mm and |2;0pv' x sin{6py)\ < 1.5 mm 
respectively 

• probability > 0.01 for reconstructed tracks with pt > 10 GeV, to remove 
high Pt miss-measured tracks 

This provides an alternative method to measure the internal structure of 
jets based on charged particles. The measurements are limited to jets with 
\y\ < 1.9, as dictated by the tracking coverage and the chosen size of the jet. 
Figures [4.231 and 14.241 show the differential jet shapes in this rapidity range using 
tracks and clusters respectively. The ratio of the jet shapes using clusters and 
tracks in data are compared to that in the Monte Carlos simulation (Figure 14. 25p . 
Maximum deviations from one of these ratios are about 5%, well within the 
quoted systematic uncertainties. 

Finally, the jet shapes results are obtained using calorimeter towers of fixed 
size 0.1 X 0.1 {y — <t> space) instead of topological clusters as input to the jet 
reconstruction algorithm. In order to suppress noise, the towers are built only 
with cells that belong to a cluster. Figure 14.261 show that tower jets are signif- 
icantly broader than cluster jets before correcting for detector effects. Instead, 
differences reduce to less than 5% after applying the corrections, as shown in Fig- 
ure 14.271 Again, the observed differences are well inside the quoted uncertainty 
on the nominal results based on calorimeter clusters. 
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Figure 4.23: Measured differential jet shapes using tracks inside jets for jets with 
\y\ < 1.9 and 30 GeV < pr < 210 GeV. 
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Figure 4.24: Measured differential jet siiapes using calorimeter clusters for jets with 
\y\ < 1.9 and 30 GeV < pr < 210 GeV. 
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Figure 4.25: Double-ratio of the differential jet shapes derived by comparing calorime- 
ter and tracking ratios of results in data and Monte Carlo simulations for jets with 
\y\ < 1.9 and 30 GeV < pr < 210 GeV. 
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Figure 4.26: Comparison of differential jet sliapes witli calorimeter clusters and topo- 
towers before correcting for detector effects for jets with \y\ < 1.9 and 30 GeV < pT < 
210 GeV. 
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Figure 4.27: Comparison of differential jet shapes witli calorimeter clusters and topo- 
towers after correcting for detector effects for jets with \y\ < 1.9 and 30 GeV < pT < 
210 GeV. 
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4.7 Results 

The measurements presented in this Chapter refer to differential and integrated 
jet shapes, p(r) and \l/(r), corrected at the particle level and obtained for anti-kt 
jets with distance parameter R = 0.6 in the region \y\ < 2.8 and 30 GeV < 
Pt < 600 GeV. The measurements are presented in separate bins of pt and \y\. 
Tabulated values of the results are available in Ref. [65] . 

Figures l4.28l to fl4.31l to 14. 33p show the measured differential (integrated) 
jet shapes as a function of r in different px ranges. The dominant peak at small 
r indicates that the majority of the jet momentum is concentrated close to the 
jet axis. At low pt, more than 80% of the transverse momentum is contained 
within a cone of radius r = 0.3 around the jet direction. This fraction increases 
up to 95% at very high p^, showing that jets become narrower as pt increases. 
This is also observed in Figure H.34[ where the measured 1 — \l/(0.3), the fraction 
of the jet transverse momentum outside a fixed radius r = 0.3, decreases as a 
function of px- 

The data are compared to predictions from HERWIG++, ALPGEN, PYTHIA- 
Perugia2010, and PYTH1A-MC09 in Figure S^Hl to Fig.^Mia); and to predic- 
tions from PYTHIA-DW and PYTH1A-Perugia2010 with and without UE con- 
tributions in Figure H.34r b). The jet shapes predicted by PYTHIA-Perugia2010 
provide a reasonable description of the data, while HERWIG++ predicts broader 
jets than the data at low and very high pt- The PYTHIA-DW predictions 
are in between PYTHIA-Perugia2010 and HERWIG++ at low pt and produce 
jets which are slightly narrower at high pt- ALPGEN is similar to PYTHIA- 
Perugia2010 at low pr, but produces jets significantly narrower than the data at 
high Pt- PYTHIA-MC09 tends to produce narrower jets than the data in the 
whole kinematic range under study. The latter may be attributed to an inad- 
equate modeling of the soft gluon radiation and UE contributions in PYTHIA- 
MC09 samples, in agreement with previous observations of the particle flow 
activity in the final state [66]. Finally, Figure I4.34r b) shows that PYTHIA- 
Perugia2010 without UE contributions predicts jets much narrower than the data 
at low Pt- This confirms the sensitivity of jet shape observables in the region 
Pt < 160 GeV to a proper description of the UE activity in the final state q 

The dependence on \y\ is shown in Figure 351 and r4.36[ where the measured 
jet shapes are presented separately in five different jet rapidity regions and dif- 

■^These conclusions are supported by energy flow results at detector level, reported in Ap- 
pendix [Cj 
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ferent pt bins, for jets with px < 500 GeV. At high px, the measured 1 — \[^(0.3) 
shape presents a mild \y\ dependence, indicating that the jets become slightly 
narrower in the forward regions. This tendency is observed also in the various 
MC samples. Similarly, Figures 14.371 and [4. 38l present the measured 1 — \E'(0.3) as 
a function of px in the different \y\ regions compared to PYTHIA-Perugia2010 
predictions. 

Finally, and only for illustration, the typical shapes of quark- and gluon- 
initiated jets, as determined using events generated with PYTHIA-Perugia2010, 
are also shown in Figures l4.37l and 14.381 For this purpose, MC events are selected 
with at least two particle-level jets with p^ > 30 GeV and \y\ < 2.8 in the final 
state. The two leading jets in this dijet sample are classified as quark-initiated or 
gluon-initiated jets by matching {in y — (p space) their direction with one of the 
outgoing partons from the QCD 2 — ?■ 2 hard process. At low p^ the measured 
jet shapes are similar to those from gluon-initiated jets, as expected from the 
dominance of hard processes with gluons in the final state. At high pt, where 
the impact of the UE contributions becomes smaller (see Figure I4.34r b)). the 
observed trend with pt in the data is mainly attributed to a changing quark- and 
gluon-jet mixture in the final state, convoluted with perturbative QCD effects 
related to the running of the strong coupling. 

The jet shapes results presented in this Chapter were obtained with the first 
pp collisions at -y/i = 7 TeV, in which the pile-up can be neglected, and were 
published in [59]. They indicate the potential of jet shape measurements at 
the LHC to constrain the phenomenological models for soft gluon radiation, 
underlying event (UE) activity, and non-perturbative fragmentation processes in 
the final state. After the publication of these results, additional comparisons 
to predictions from different MC programs have been performed [67]. They are 
reported in the next Chapter, which includes tests to the data points for all 
the MC models considered. 
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Figure 4.28: The measured differential jet shape, p{r), in inclusive jet production for 
jets with |y| < 2.8 and 30 GeV < pT < 110 GeV is shown in different px regions. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. 
The predictions of PYTHIA-Perugia2010 (sohd lines), HERWIG++ (dashed hues), 
ALPGEN interfaced with HERWIG and JIMMY (dotted hues), and PYTHIA-MC09 
(dashed-dotted lines) are shown for comparison. 
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Figure 4.29: The measured differential jet shape, p{r), in inclusive jet production for 
jets with \y\ < 2.8 and 110 GeV < pT < 310 GeV is shown in different px regions. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. 
The predictions of PYTHIA-Perugia2010 (sohd lines), HERWIG++ (dashed hues), 
ALPGEN interfaced with HERWIG and JIMMY (dotted hues), and PYTHIA-MC09 
(dashed-dotted lines) are shown for comparison. 
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Figure 4.30: The measured differential jet shape, p{r), in inclusive jet production for 
jets with \y\ < 2.8 and 310 GeV < pT < 600 GeV is shown in different px regions. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. 
The predictions of PYTHIA-Perugia2010 (sohd lines), HERWIG++ (dashed hues), 
ALPGEN interfaced with HERWIG and JIMMY (dotted hues), and PYTH1A-MC09 
(dashed-dotted lines) are shown for comparison. 
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Figure 4.31: The measured integrated jet shape, ^'(r), in inclusive jet production for 
jets with \y\ < 2.8 and 30 GeV < pT < 110 GeV is shown in different px regions. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. 
The predictions of PYTHIA-Perugia2010 (sohd lines), HERWIG++ (dashed hues), 
ALPGEN interfaced with HERWIG and JIMMY (dotted hues), and PYTHIA-MC09 
(dashed-dotted lines) are shown for comparison. 
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Figure 4.32: The measured integrated jet shape, ^(r), in inclusive jet production for 
jets with \y\ < 2.8 and 110 GeV < pT < 310 GeV is shown in different px regions. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. 
The predictions of PYTHIA-Perugia2010 (sohd lines), HERWIG++ (dashed hues), 
ALPGEN interfaced with HERWIG and JIMMY (dotted hues), and PYTHIA-MC09 
(dashed-dotted lines) are shown for comparison. 
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Figure 4.35: The measured integrated jet shape, 1 — ^'(r = 0.3), as a function 
of \y\ for jets with \y\ < 2.8 and 30 GeV < px < 210 GeV. Error bars indicate 
the statistical and systematic uncertainties added in quadrature. The predictions of 
PYTHIA-Perugia2010 (solid lines), HERWIG++ (dashed lines), ALPGEN interfaced 
with HERWIG and JIMMY (dotted hues), PYTHIA-MC09 (dashed-dotted lines), and 
PYTHIA-DW (daslied-dotted-dotted lines) are shown for comparison. 
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Figure 4.36: The measured integrated jet shape, 1 — ^(r = 0.3), as a function of 
\y\ for jets with \y\ < 2.8 and 210 GeV < px < 500 GeV. Error bars indicate 
the statistical and systematic uncertainties added in quadrature. The predictions of 
PYTHIA-Perugia2010 (soHd lines), HERWIG++ (dashed lines), ALPGEN interfaced 
with HERWIG and JIMMY (dotted hues), PYTHIA-MC09 (dashed-dotted lines), and 
PYTHIA-DW (dashed-dotted-dotted lines) are shown for comparison. 
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Figure 4.37: The measured integrated jet shape, 1 — ^{r = 0.3), as a function of px 
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Figure 4.38: The measured integrated jet shape, 1 — ^{r = 0.3), as a function oi px in 
different jet rapidity regions for jets with \y\ < 2.8 and 30 GeV < px < 500 GeV. Error 
bars indicate the statistical and systematic uncertainties added in quadrature. The 
predictions of PYTHIA-Perugia2010 (solid line) are shown for comparison, together 
with the prediction separately for quark-initiated (dashed lines) and gluon-initiated 
jets (dotted lines) in dijet events. 
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Chapter 5 



Jet Shapes in ATLAS and Monte 
Carlo modeling 

In this Chapter the jet shapes results pubhshed in [59] are compared to predic- 
tions of different MC programs, that include recent tuned sets of parameters to 
describe the ATLAS minimum bias (MB) and jet data. 

5.1 Monte Carlo samples 

The MC generators considered in this Chapter are listed below. 

• The PYTHIA 6.4.23 generator is used to produce QCD-jet inclusive sam- 
ples with both AMBTl and AUET2 [SS] tunes, as determined by ATLAS, 
and with the Perugia2011 tune [20]. PYTHIA-AMBTl is based on MRST 
LO* parton density functions (PDFs) inside the proton, and uses as input 
results from the ATLAS charged particle MB analysis at ^/s = 0.9 TeV 
and = 7 TeV. PYTHIA-AUET2 is based on MRST LO** PDFs, uses 
also as input ATLAS UE data [69], and includes a modified final-state par- 
ton shower from initial-state radiated partons to describe the CDF [53j 
and ATLAS jet shape measurements. The PYTHIA- Perugia2011 tune 
uses CTEQ5L PDFs and constitutes an updated version of the PYTHIA- 
Perugia2010 tune, as inspired by the comparison with early LHC data. 
Samples are also generated using the PYTHIA 8.145 program [70] which, 
among other updates, includes a fully interleaved px-ordered evolution for 
initial- and final-state gluon radiation and multiple parton interactions 
(MPI). The Pythia 8 samples are produced with the 4C tune [TD] for the 
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UE modeling, based on CTEQ6L1 PDFs and using Tevatron and ATLAS 
data as input. 

• QCD-jet inclusive samples are generated with the Herwig++ 2.4.2 and 
Herwig++ 2.5.1 programs. The Herwig++ 2.4.2 predictions used in [59] 
suffered from a wrong setting in the ATLAS generation of the MC sam- 
ples, related to the parameters that control the contribution from MPI, 
which a priori would affect the comparison with the data at low pt- In 
this Chapter, the same Herwig++ 2.4.2 predictions already presented in 
the previous Chapter (here denoted as Herwig++ 2.4.2 bug), are compared 
with correctly generated HERWIG 2.4.2 predictions. Herwig++ 2.5.1 in- 
cludes a modified algorithm for color reconnection between clusters before 
hadronization and new color connections between soft scatters and the 
beam remnants that improve the description of the MB data and UE- 
related observables at the LHC. In the case of Herwig++ 2.4.2 the LO* 
tune [7T] is employed, while the UE7 tune [72] is used for Herwig++2.5.1 
predictions. In addition, samples are produced with HERWIG 6.510 inter- 
faced with JIMMY 4.31 to model the UE contributions with both AUETl- 
LO* [73] and AUET2-L0** [68] tunes, as determined by ATLAS using data 
at 0.9 TeV and 7 TeV. 

• QCD-jet events are generated using the ALPGEN 2.13 program interfaced 
with PYTHIA 6.423 for the parton shower and fragmentation into hadrons. 
The samples are produced with CTEQ6L1 PDFs and the Perugia6 tune to 
model the UE contributions. 

• The Sherpa 1.2.3 and Sherpa 1.3.0 programs [73] are used to generate 
QCD-jet samples. Default settings include CTEQ66 PDFs and MPI tuned 
to Tevatron and ATLAS data. The Sherpa 1.3.0 version contains a modi- 
fied parton shower evolution via pre-factors within as{pT), as motivated by 
next-to-leading-logarithm (NLL) perturbative QCD (pQCD) calculations, 
that improve the description of the measured Z boson pt in ATLAS. Sam- 
ples with different matrix element multiplicities are produced, with the aim 
to illustrate the impact of higher-order contributions to the description of 
the jet shapes. Samples with only 2 — )■ 2 hard processes and with up to 
2 — 7- 6 partons in the final state are considered, and interfaced with parton 
showers. 

• QCD-jet samples are generated with the POWHEG-BOX r302 program. 
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which includes next-to-leading order (NLO) pQCD predictions for inclu- 
sive jet production, interfaced either with PYTHIA 6.423 or with HER- 
WIG 6.510 plus JIMMY 4.31 for parton shower, fragmentation into hadrons, 
and to model the underlying event. For the latter, the AMBTl or AUETl- 
LO* tunes are used, respectively. 

The simulated event samples were produced using the ATLAS interface to the 
generators and the Rivet system was used to analyze the samples and produce 
the relevant jet shape observables. 

5.2 Results 

Different observables in data are compared to the MC predictions. This includes 
the measured differential jet shape p{r) for jets with \y\ < 2.8 in eight pt bins 
within the range 30 < pt < 310 GeV, and the integrated jet shape 1 — \E'(0.3) 
as a function of \y\ in different pt regions and as a function of pt in different \y\ 
bins. 

5.2.1 Comparison with PYTHIA 

The measured jet shapes are compared to the different PYTHIA predictions 
in Figures 15.11 to 15.51 As expected, PYTHIA6- AMBTl tends to produce jets 
slightly narrower than the data as it underestimates the UE activity in dijet 
events [68] and lacks a tuned final-state parton shower from the initial-state 
radiation. PYTHIA6-AUET2, PYTHIA6-Perugia2011, and Pythia 8-4C provide 
a good description of the measured jet shapes in the whole jet kinematic range 
considered 0. 

5.2.2 Comparison with Herwig++ and HERWIG /JIMMY 

The measurements are compared to the predictions from Herwig++ 2.4.2, Her- 
wig++ 2.5.1, and HERWIG/JIMMY programs in Figure ESI to Figure IHTTOl As 
shown in the Figures, the effect of the wrong setting in Herwig++ 2.4.2 bug (see 
Section 2) is only visible for jets with px < 40 GeV. Herwig++ 2.4.2 provides a 
reasonable description of the data in the whole jet kinematic range considered. 

^The reader should note that the ATLAS jet shapes measurements are input to the 
PYTHIA6-AUET2 tune. 
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HERWIG/ JIMMY, which includes the AUET2 tune, describes the data best. 
Finally, Herwig++ 2.5.1 tends to produce jets narrower than the data. 

5.2.3 Comparison with ALPGEN and Sherpa 

The data are compared to the predictions from Sherpa 1.2.3 (2 — )■ 2 process), 
Sherpa 1.2.3 (up to 2 — )■ 6 process), Sherpa 1.3.0 (2 — )■ 2 process), and ALPGEN 
interfaced with PYTHIA in Figure 15.111 to Figure 15.151 Comparisons with ALP- 
GEN interfaced with HERWIG and JIMMY were presented in [59]. The different 
Sherpa predictions are similar and provide a reasonable description of the data. 
This indicates that the presence of additional partons from higher-order matrix 
elements contributions do not affect the predicted jet shapes, mainly dictated by 
the soft radiation in the parton shower. The comparison between Sherpa 1.2.3 
and Sherpa 1.3.0 shows that the NLL-inspired corrections included in the latter 
for the parton shower do not impact significantly the predicted jet shapes. ALP- 
GEN interfaced with PYTHIA predicts too-narrow jets and does not describe 
the data. This was already the case for ALPGEN interfaced with HERWIG 
and JIMMY [59] and requires further investigation to determine whether the 
disagreement observed with the data can be completely attributed to the UE 
modeling in the MC samples or it is also related to the prescription followed 
by ALPGEN in merging the partons from the matrix elements with the parton 
showers in the final state. 

5.2.4 Comparison with POWHEG 

The impact of the presence of higher-order matrix elements contributions on the 
predicted jet shapes is further studied in Figure 15.161 to Figure I5.20[ where the 
data are compared to the predictions from POWHEG interfaced with PYTHIA 
or HERWIG/ JIMMY, as well as to standalone PYTHIA and HERWIG/ JIMMY 
predictions. For these comparisons, AMBTl and AUETl-LO* tuned sets of pa- 
rameters are used for PYTHIA and HERWIG/ JIMMY, respectively. In general, 
the jet shapes predicted by POWHEG follow the trend of the corresponding 
standalone predictions. POWHEG interfaced with HERWIG/ JIMMY provides 
a reasonable description of the data while the interface with PYTHIA predicts 
too-narrow jets. This confirms that the shape of the jet is mainly dictated by the 
parton shower implementation and the details of the UE modeling in the final 
state. 
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Figure 5.1: The measured differential jet shape, p{r), in inclusive jet production for 
jets with \y\ < 2.8 and 30 GeV < px < 110 GeV is shown in different pT regions. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. 
The predictions of PYTHIA-Perugia2011 (sohd lines), PYTHIA-AUET2 (dotted lines), 
PYTHIA-AMBTl (dashed lines), and Pythia 8-4C (dashed-dotted lines) are shown for 
comparison. 
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Figure 5.2: The measured differential jet shape, p{r), in inclusive jet production for 
jets with \y\ < 2.8 and 110 GeV < pr < 310 GeV is shown in different px regions. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. 
The predictions of PYTHIA-Perugia2011 (sohd lines), PYTHIA-AUET2 (dotted lines), 
PYTHIA-AMBTl (dashed lines), and Pythia 8-4C (dashed-dotted lines) are shown for 
comparison. 
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Figure 5.3: The measured integrated jet shape, 1 — ^{r = 0.3), as a function of 
\y\ in different jet px regions for jets with \y\ < 2.8 and 30 GeV < px < 110 GeV. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. 
The predictions of PYTHIA-Perugia2011 (sohd hues), PYTHIA- AUET2 (dotted hues), 
PYTHIA- AMBTl (dashed lines), and Pythia 8-4C (dashed-dotted lines) are shown for 
comparison. 
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Figure 5.4: The measured integrated jet shape, 1 — ^'(r = 0.3), as a function of \y\ 
in different jet pT regions for jets with \y\ < 2.8 and 110 GeV < px < 310 GeV. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. 
The predictions of PYTHIA-Perugia2011 (sohd Unes), PYTHIA- AUET2 (dotted hues), 
PYTHIA- AMBTl (dashed lines), and Pythia 8-4C (dashed-dotted hues) are shown for 
comparison. 
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Figure 5.5: The measured integrated jet shape, 1 — ^{r = 0.3), as a function of px 
in different jet rapidity regions for jets with \y\ < 2.8 and 30 GeV < px < 500 GeV. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. 
The predictions of PYTHIA-Perugia2011 (sohd hues), PYTHIA- AUET2 (dotted hues), 
PYTHIA-AMBTl (dashed lines), and Pythia 8-4C (dashed-dotted hues) are shown for 
comparison. 
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Figure 5.6: The measured differential jet shape, p{r), in inclusive jet production for 
jets with \y\ < 2.8 and 30 GeV < px < 110 GeV is shown in different pT regions. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. 
The predictions of Herwig++2.4.2 (solid lines), Herwig++2.4.2 bug (dotted lines), 
Herwig++ 2.5.1 (dashed-dotted lines), and HERWIG/ JIMMY- AUET2 (dashed lines) 
are shown for comparison. 
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Figure 5.7: The measured differential jet shape, p{r), in inclusive jet production for 
jets with \y\ < 2.8 and 110 GeV < pr < 310 GeV is shown in different pT regions. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. 
The predictions of Herwig++2.4.2 (solid lines), Herwig++2.4.2 bug (dotted lines), 
Herwig++ 2.5.1 (dashed-dotted lines), and HERWIG/ JIMMY- AUET2 (dashed lines) 
are shown for comparison. 
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Figure 5.8: The measured integrated jet shape, 1 — ^{r = 0.3), as a function of 
\y\ in different jet px regions for jets with \y\ < 2.8 and 30 GeV < px < 110 GeV. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. 
The predictions of Herwig++2.4.2 (solid lines), Herwig++2.4.2 bug (dotted hues), 
Herwig++ 2.5.1 (dashed-dotted hues), and HERWIG/ JIMMY- AUET2 (dashed hues) 
are shown for comparison. 
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Figure 5.9: The measured integrated jet shape, 1 — ^'(r = 0.3), as a function of \y\ 
in different jet pT regions for jets with \y\ < 2.8 and 110 GeV < px < 310 GeV. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. 
The predictions of Herwig++2.4.2 (sohd hues), Herwig++2.4.2 bug (dotted hues), 
Herwig++ 2.5.1 (dashed-dotted hues), and HERWIG/ JIMMY- AUET2 (dashed hues) 
are shown for comparison. 
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Figure 5.10: The measured integrated jet shape, 1 — ^{r = 0.3), as a function of px 
in different jet rapidity regions for jets with \y\ < 2.8 and 30 GeV < pT < 500 GeV. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. 
The predictions of Herwig++2.4.2 (sohd hues), Herwig++2.4.2 bug (dotted hues), 
Herwig++ 2.5.1 (dashed-dotted hues), and HERWIG/ JIMMY- AUET2 (dashed hues) 
are shown for comparison. 
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Figure 5.11: The measured differential jet shape, p{r), in inclusive jet production for 
jets with \y\ < 2.8 and 30 GeV < px < 110 GeV is shown in different pT regions. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. 
The predictions of Sherpa 1.3.0 (2 — )• 2)(solid lines), Sherpa 1.2.3 (2 — )• 2) (dotted 
lines), Sherpa (up to 2 — )• 6) (dashed lines), and ALPGEN interfaced to PYTHIA 
(dashed-dotted lines) are shown for comparison. 
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Figure 5.12: The measured differential jet shape, p{r), in inclusive jet production for 
jets with \y\ < 2.8 and 110 GeV < pT < 310 GeV is shown in different pj- regions. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. 
The predictions of Sherpa 1.3.0 (2 — )• 2)(solid lines), Sherpa 1.2.3 (2 — )• 2) (dotted 
lines), Sherpa (up to 2 — >• 6) (dashed lines), and ALPGEN interfaced to PYTHIA 
(dashed-dotted lines) are shown for comparison. 
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Figure 5.13: The measured integrated jet shape, 1 — ^{r = 0.3), as a function of 
\y\ in different jet px regions for jets with \y\ < 2.8 and 30 GeV < px < 110 GeV. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. 
The predictions of Sherpa 1.3.0 (2 — > 2)(sohd hues), Sherpa 1.2.3 (2 — )• 2) (dotted 
hues), Sherpa (up to 2 — )• 6) (dashed hnes), and ALPGEN interfaced to PYTHIA 
(dashed-dotted hnes) are shown for comparison. 
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Figure 5.14: The measured integrated jet shape, 1 — ^{r = 0.3), as a function of 
\y\ in different jet px regions for jets with \y\ < 2.8 and 110 GeV < pT < 310 GeV. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. 
The predictions of Sherpa 1.3.0 (2 — > 2)(solid hues), Sherpa 1.2.3 (2 — )• 2) (dotted 
lines), Sherpa (up to 2 — )• 6) (dashed lines), and ALPGEN interfaced to PYTHIA 
(dashed-dotted lines) are shown for comparison. 
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Figure 5.15: The measured integrated jet shape, 1 — "^{r = 0.3), as a function of pT 
in different jet rapidity regions for jets with \y\ < 2.8 and 30 GeV < pT < 500 GeV. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. 
The predictions of Sherpa 1.3.0 (2 — )• 2)(sohd hues), Sherpa 1.2.3 (2 — )• 2) (dotted 
hues), Sherpa (up to 2 — )• 6) (dashed hues), and ALPGEN interfaced to PYTHIA 
(dashed-dotted hues) are shown for comparison. 
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Figure 5.16: The measured differential jet siiape, p{r), in inclusive jet production for 
jets with \y\ < 2.8 and 30 GeV < pr < 110 GeV is shown in different px regions. Er- 
ror bars indicate the statistical and systematic uncertainties added in quadrature. The 
predictions of POWHEG interfaced with PYTHIA-AMBTl (dashed hues), POWHEG 
interfaced with HERWIG/ JIMMY- AUETl (soUd Hues), PYTHIA-AMBTl (dotted 
lines), and HERWIG/ JIMMY- AUETl (dashed-dotted lines) are shown for compar- 



ison. 
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Figure 5.17: The measured differential jet siiape, p{r), in inclusive jet production for 
jets with \y\ < 2.8 and 110 GeV < pT < 310 GeV is shown in different pT regions. Er- 
ror bars indicate the statistical and systematic uncertainties added in quadrature. The 
predictions of POWHEG interfaced with PYTHIA-AMBTl (dashed hues), POWHEG 
interfaced with HERWIG/ JIMMY- AUETl (sohd hues), PYTHIA-AMBTl (dotted 
lines), and HERWIG/ JIMMY- AUETl (dashed-dotted lines) are shown for compar- 



ison. 
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Figure 5.18: The measured integrated jet shape, 1 — ^'(r = 0.3), as a function of \y\ in 
different jet pT regions for jets with \y\ < 2.8 and 30 GeV < pT < 110 GeV. Error bars 
indicate the statistical and systematic uncertainties added in quadrature. The predic- 
tions of POWHEG interfaced with PYTHIA-AMBTl (dashed hues), POWHEG inter- 
faced with HERWIG/JIMMY-AUETl (sohd hues), PYTHIA-AMBTl (dotted hues), 
and HERWIG/JIMMY-AUETl (dashed-dotted hues) are shown for comparison. 
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Figure 5.19: The measured integrated jet shape, 1 — ^'(r = 0.3), as a function of \y\ in 
different jet px regions for jets with \y\ < 2.8 and 110 GeV < pT < 310 GeV. Error bars 
indicate the statistical and systematic uncertainties added in quadrature. The predic- 
tions of POWHEG interfaced with PYTHIA-AMBTl (dashed hues), POWHEG inter- 
faced with HERWIG/JIMMY-AUETl (sohd hues), PYTHIA-AMBTl (dotted hues), 
and HERWIG/JIMMY-AUETl (dashed-dotted hnes) are shown for comparison. 
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Figure 5.20: The measured integrated jet shape, 1 — ^'(r = 0.3), as a function of in 
different jet rapidity regions for jets with \y\ < 2.8 and 30 GeV < px < 500 GeV. Er- 
ror bars indicate the statistical and systematic uncertainties added in quadrature. The 
predictions of POWHEG interfaced with PYTHIA-AMBTl (dashed hues), POWHEG 
interfaced with HERWIG/ JIMMY- AUETl (sohd hues), PYTHIA-AMBTl (dotted 
lines), and HERWIG/ JIMMY- AUETl (dashed-dotted lines) are shown for compar- 
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Figure 5.21: The measured integrated jet shape, 1 — ^{r = 0.3), as a function of px 
for jets with \y\ < 2.8 and 30 GeV < px < 600 GeV. Error bars indicate the statistical 
and systematic uncertainties added in quadrature. The measurements are compared 
to the different MC predictions considered. 
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5.2.5 statistical tests 

Finally, a test is performed to the data points in Figures I4.38[ 15. 5[ I5.10[ I5.15[ 
and l5.20l with respect to a given MC prediction, separately in each rapidity region. 
The systematic uncertainties are considered independent and fully correlated 
across pr bins, and the test is carried out according to the formula 

= ' fjt~7/^'^'lL + h^^]" ' (5-1) 

where dj is the measured data point in bin j, mcj{s) is the corresponding MC 
prediction, and s denotes the vector of standard deviations, Si, for the different 
independent sources of systematic uncertainty. For each rapidity region consid- 
ered, the sums above run over the total number of data points in and five 
independent sources of systematic uncertainty, and the is minimized with re- 
spect to s. Correlations among systematic uncertainties are taken into account 
in mcj{s). 

The results for the different MC predictions shown in this and in the 
previous Chapter are collected in Table 5.1. 
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Table 5.1: Results of tests to the data in Figures I4.38[ 15.51 15.10[ I5.15[ and 
15.201 with respect to the different MC predictions. 
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Conclusions 



This Thesis describes the measurements of the inclusive jet cross section and jet 
shapes in pp coUisions at = 7 TeV using the ATLAS detector. The anti-fcf 
algorithm with i? = 0.6 is used to reconstruct jets. Both measurements are un- 
folded back to the hadron level and compared to theory predictions. 

The inclusive jet cross section measurement has been performed using jets with 
Pt > 20 GeV and \y\ < 4.4 in a data sample corresponding to a total luminos- 
ity of 37 pb~^ collected by ATLAS during 2010. The measurement constitute a 
stringent test of pQCD over ten orders of magnitude, covering a jet px range up 
to 1.5 TeV. The measurement expands by a factor of ~ 2 the kinematic region 
in Pt probed by the Tevatron, in a wider y range. 

The data are compared with NLO predictions using several PDFs corrected for 
non-perturbative effects. A general agreement is found, but the predicted cross 
section tends to be larger than that in data at high px in the forward region. 
This can be atributted to the current PDF of the gluon at high x, showing the 
potential of the measurement to constrain the parton distribution functions when 
used in global fits. 

The measurement of jet shapes in inclusive jet production uses 3 pb~^ collected 
by ATLAS, for jets with pr > 30 GeV and \y\ < 2.8. It is important to notice 
that the measurement was performed with the very first pp collisions at 7 TeV 
center-of-mass-energy, when the LHC luminosity was low enough to collect a 
data sample clean of pile-up. The jet shape measurements can be used to con- 
strain the phenomenological models for soft gluon radiation, underlying event 
(UE) activity, and non-perturbative fragmentation processes in the final state. 
This is particularly important for searches of new physics that require a good 
control of the soft physics. The jet shape measurements constitute a first step 
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towards the study of boosted topologies corresponding to the decay of heavy 
particles, that promise to play a central role in the serach for a hght Higgs boson 
or in searches for new physics. 

The integrated luminosity recorded by ATLAS has already increased above 1 fb^^ 
during 2011. The inchisivc jet cross section is being updated with these new data, 
and further studies are being done in order to reduce the JES uncertainty. The 
large dataset being recorded will also allow to perform sensitive jet substructure 
analyses searching for new physics. 
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Sensitivity of the jet shapes to 
the underlying event and to the 
pile-up using MC simulated 
events 

Before the first pp collisions were delivered by the LHC, studies on jet shapes 
were performed using PYTHIA QCD-dijet Monte Carlo simulated events in pp 
collisions at i/i =10 TeV. The effect of the underlying event and the pileup on 
the jet shapes was investigated in jets reconstructed with the SISCone algorithm 
with R = 0.7 and splitting-merging fraction equal to 0.75. 

Figures lA.H IA.2I and IA.3I show the comparison of jet shapes at hadron level 
in events with and without underlying event. As expected, the presence of the 
underlying event translates into broader jets at low pt- The jet shapes results 
obtained in ATLAS Monte Carlo simulated events are compared with published 
results from the CDF experiment using a midpoint algorithm, corrected back to 
the hadron level. Given the difference in the quark/gluon-jet mixture, the dif- 
ferent relative importance of the underlying event contributions in the Tevatron 
and the LHC energies and the differences in the jet reconstruction algorithms 
employed, one should not expect an a priori agreement between experiments and, 
therefore, no conclusions are made. However, it is remarkable to observe that 
SISCone simulated jet shapes in ATLAS are so similar to those published by 
CDF. 

The sensitivity of the jet shapes to pile-up contributions is studied using a 
sample of QCD-dijet events at ^/s = 10 TeV, overlaid with minimum bias inter- 
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actions corresponding to an instantaneous luminosity of 2 x 10'^^ cm~^s~^, for 
which on average 4 additional interactions are expected. The contribution from 
cavern background events, arising from neutrons that interact before thermali- 
sation, is also included. Figures IA.4[ IA.5I and IA.6I show the comparison of jet 
shapes at detector level in events with and without pile-up. Jets in events with 
pile-up are significantly broader, mainly at low pt- 
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Figure A.l: Differential jet siiapes for jets witli 37 GeV < pt < 148 GeV and 0.1 < 
\rj\ < 0.7 as reconstructed using tiie SISCone algorithm for events with (full dots) 
and without (open dots) UE contributions. In addition, results from CDF data (full 
triangles) are shown. 
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Hadron Level, Inclusive Jets 

ATLAS Pythia: SISCone (R = 0.7, SMF = 0.75) 
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Figure A. 2: Differential jet shapes for jets with 148 GeV < px < 380 GeV and 
0.1 < \r]\ < 0.7 as reconstructed using the SISCone algorithm for events with (full 
dots) and without (open dots) UE contributions. In addition, results from CDF data 
(full triangles) are shown. 
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Appendix A. Sensitivity of the jet shapes to the underlying event and 
to the pile-up using MC simulated events 



Hadron Level, Inclusive Jets 

ATLAS Pythia: SISCone (R = 0.7, SMF = 0.75) 
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Figure A. 3: Integrated jet shapes for jets with 37 GeV < pr < 380 GeV and 0.1 < 
\ri\ < 0.7 as reconstructed using the SISCone algorithm for events with (full dots) 
and without (open dots) UE contributions. In addition, results from CDF data (full 
triangles) are shown. 
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Detector Level, Leading Jet 
SISCone, R = 0.7, SMF = 0.75 



10-' 



■ ^ 37 GeV/c < p'" < 45 GeV/c 
;A A A A Jl^"l<l-2 

" ^ A ^ 

A A 

A Pythia w/o Pile-up ^ 
A Pythia with Hle-up 

; ... 1 .... 1 .... 1 .... 1 .... 1 .... 1 .... 1 .... 1 .... 1 ... . 


■ ^ 97 Ge V/c < p^J:' < 1 1 2 Ge V/c 

A lTf°'l<1.2 

r ^ A 

A A 

A 
A 

: ... 1 .... 1 .... 1 .... 1 .... 1 .... 1 .... 1 .... 1 .... 1 ... . 


■ * 186 GeV/c < p'" < 208 GeV/c 
4 lif'l<1.2 

A 

A A 

- 4 

: ... 1 .... 1 .... 1 .... 1 .... 1 .... 1 .... 1 .... 1 .... 1 ... . 


■ * 340 GeV/c < p'" < 380 GeV/c 

^ lTf"l<1.2 

4 

' 4 , 

■ ■ ■ ■ 1 1 .... 1 .... 1 .... 1 .... 1 .... 1 .... 1 ... . 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



r/R 



Figure A. 4: Differential jet shapes for the leading jet with 37 GcV < pT < 380 GeV 
and \r)\ < 1.2 events with (open triangles) and without (full triangles) pile-up. 
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Figure A. 5: Integrated jet shapes for the leading jet with 37 GeV < px < 380 GeV 
and \r]\ < 1.2 events with (open triangles) and without (full triangles) pile-up. 
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Figure A.6: 1 - ^'(0.3) for the leading jet with 37 GeV < pr < 380 GeV and |r/| < 1.2 
events with (open triangles) and without (full triangles) pile-up. 
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Appendix A. Sensitivity of the jet shapes to the underlying event and 
to the pile-up using MC simulated events 
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Appendix B 



Jet shapes and energy flow in pp 
collisions at y/s = 900 GeV 

In this appendix, the measurement of jet shapes and energy flow at detector 
level in pp collisions at y/s = 900 GeV is presented. The data are compared to 
different Monte Carlo predictions including full ATLAS detector simulation. 

B.l Event Selection and Monte Carlo Simula- 
tion 

The analysis is based on a minimum bias (MB) data sample collected by the 
ATLAS experiment in 2009. Data was selected from luminosity blocks where 
the inner detector (ID) was in fully operational state, the data quality (DQ) 
flags for the calorimeter were good, and only events in filled and paired bunch 
counter BCIDs were accepted. The Minimum Bias trigger was employed. Events 
are then required to have a reconstructed primary vertex with 2;-position within 
10 cm of the nominal interaction point and with at least three tracks pointing to 
it. Jets are reconstructed from the energy deposits in the calorimeter using the 
anti-/ct algorithm with D = 0.6 and topotowers as input. 

The events are required to have at least a jet with uncorrected transverse 
momentum p^ above 7 GeV and rapidity in the range \y\ < 2.6. In addition, 
jets affected by noise in the calorimeter are rejected. The measurements are 
compared to Monte Carlo MB simulated events as generated using pythia and 
PHOJET [75] event generators, interfaced with a full simulation of the ATLAS 
detector response to particles based on geant4. In the case of pythia, different 
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Appendix B. Jet shapes and energy flow in pp collisions at ^/s = 900 GeV 



Monte Carlo samples with slightly different tunes for the parton shower and 
underlying event modeling in the final state have been considered. Samples are 
generated with ATLAS MC09, DW and PerugiaO tunes. The simulated events 
are passed through the same trigger selection and analysis chain as in the data. 

Figure IB. II shows the measured jet multiplicity, px, y and distributions 
compared to predictions from MB Monte Carlo events (pythia tune ATLAS 
MC09). The measurements are reasonably well described by the simulation. 
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Figure B.l: Measured basic jet kinematic distributions compared to Monte Carlo 
predictions (pythia tune ATLAS MC09), normalized to the number of jets ob- 
served in data. 



B.1.1 Jet Shapes using Calorimeter Towers 

Figure IB. 21 shows the measured differential and integrated jets shapes, as de- 
termined using calorimeter towers, for jets with pt > 7 GeV and \y\ < 2.6. 
Similarly, Figs. IB. 3 1 to [K4l presents the measurements in separate bins of jet ra- 
pidity. The profiles present the expected shape with a prominent peak at low r 
which indicates that the majority of the jet momentum is concentrated in the 
core of the jet. The measured jet shapes do not present a significant rapidity 
dependence. This is clearly observed in Fig. IB. 5 1 where the measured integrated 
jet shape \E'(r = 0.3) is presented as a function of \y\. The data are compared 
with Monte Carlo simulations from the different MB models. The predictions 
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B.1.2 Jet Shapes using Tracks 



from PYTHIA MC09, pythia PerugiaO and phojet are similar and close to the 
measurements, although predict jets slightly narrower than the data, pythia 
DW produces jets slightly broader than the data. 

B.1.2 Jet Shapes using Tracks 

The tracking system provides an alternative approach to measure the internal 
structure of the jets based on charged particles. Tracks are selected according to 
the following criteria: 

• Pt> 500 MeV and |r/| < 2.5, 

• number of associated Pixel hits > 1 

• number of associates SCT hits > 6, 

• |(iOpy| < 1.5 mm, 

• \zQpvsin{6pv)\ < 1.5 mm, 

where dOpy and zQpv are the impact parameter and z position at the perigee 
measured with respect to the primary vertex. 

Figure IB. 61 shows the track multiplicity inside the jets compared to Monte 
Carlo predictions before and after applying final track quality cuts. In both 
cases, the simulation provides a reasonable description of the data. Figure IB. 71 
shows the measured differential jet shapes and jet profiles, as determined using 
tracks, for jets with pt > 7 GeV and \y\ < 1.9. The measurements are presented 
in different rapidity bins in Fig. IB. 81 The conclusions are similar to those on the 
calorimeter quantities. The simulation provides a reasonable description of the 
data although tends to produce jets slightly narrower than the data. 
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Appendix B. Jet shapes and energy flow in pp collisions at ^/s = 900 GeV 
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Figure B.2: Measured differential and integrated jet shapes using calorimeter 
towers for jets with pt > 7 GeV and \y\ < 2.6. The data are compared to 
various Monte Carlo simulations. 
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Figure B.3: Measured differential jet shapes using calorimeter towers for jets 
with Pt > 7 GeV as a function of \y\. The data are compared to various Monte 
Carlo simulations. 
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B.1.2 Jet Shapes using Tracks 
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Figure B.4: Measured integrated jet shapes using calorimeter towers for jets with 
Pt > 7 GeV as a function of \y\. The data are compared to various Monte Carlo 
simulations. 
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Appendix B. Jet shapes and energy flow in pp collisions at ^/s = 900 GeV 
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Figure B.5: Measured integrated jet shapes \l/(r = 0.3) using calorimeter towers 
for jets with pr > 7 GeV as a function of \y\. The data are compared to various 
Monte Carlo simulations. 




Figure B.6: Measured total number of tracks inside the jet before and after 
final track quality cuts, for jets with p^ > 7 GeV and \y\ < 1.9. The data are 
compared to various Monte Carlo simulations. 



144 



B.1.2 Jet Shapes using Tracks 
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Figure B.7: Measured differential jet shapes using trades for jets with px > 7 GeV 
and \y\ < 1.9. The data are compared to various Monte Carlo simulations. 
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Figure B.8: Measured differential jet shapes using tracks for jets with px > 7 GeV 
as a function of \y\. The data are compared to various Monte Carlo simulations. 
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Appendix B. Jet shapes and energy flow in pp collisions at ^/s = 900 GeV 



B.2 Energy Flow 

The study of the particle flow away from the jet direction provides additional 
information on the soft radiation and the underlying event contribution in the 
flnal state, as well as testing the accuracy of the Monte Carlo simulation of the 
detector response to low energy particles. 

B.2.1 Energy flow in the azimuthal direction 

In this analysis, the hadronic activity outside the jet cone is studied using either 
calorimeter towers or tracks, and expressed in terms of the average transverse 
momentum observed as a function of the distance to the jet axis in the azimuthal 
direction: 

\d(j)\dy 2D\A(j)\ Njet 

(B.l) 

where — A0/2|, |0 + A0/2|) is the scalar sum of the transverse momentum 
of the constituents (calorimeter or tracks) at a given distance to the jet. Bins 
of A0 = 0.2 have been used. As illustrated in Fig lB.Qt only constituents within 
the rapidity range occupied by the jet cone are considered, with the aim to limit 
the effect of a different calorimeter response as a function of rapidity. 
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Figure B.9: Sketch of the particle flow as a function of the distance in azimuth 
to the jet axis. 

Figures IB. 101 and IB. Ill show the measured energy flow using calorimeter towers 
for events with at least one jet with px >7 GeV, in different jet rapidity regions 
up to \y\ = 2.6. The measurements are compared to the predictions from differ- 
ent MB simulated samples. The Monte Carlo simulation provides a reasonable 



146 



B.2.1 Energy flow in the azimuthal direction 



description of the core of the jet (| A(/)| < 0.6) but tends to shghtly underestimate 
the region between jets with \A(f)\ ~ 0/2, dominated by soft hadronic activity. 
As A(j) increases, the measured energy flow increases due to the presence of a sec- 
ond jet. Similarly, Figs. IB. 121 and IB. 131 show the measured flow using tracks, in 
events with jets with px >7 GeV and \y\ < 1.9, compared to various Monte Carlo 
predictions. The different pythia simulations slightly overestimate the particle 
flow at the core of the jet while slightly underestimate the activity in the region 
orthogonal to the jet direction. The latter effect is particularly significant in the 
case of PHOJET samples in both track- and calorimeter-based measurements. 



[GeV/bin] 


: ' Ant'i-k-^l b = 0.6 

- p'^" >7GeV and |y''"|<1.9 

■*H 1 yi=l _ ylower 1 ^ □ 


• 900 GeV MinBias DATA - 

900 GeV Pythia IVIC09 

900 GeV Pythia Perugia 

■■■■■900 GeV Pythia DW 
-■■■-900GeVPhojet 


(/} 




Beam 2009 






UNCORRECTED - 


-©- 
<1 




Statistical errors only : 


Q 

<M 

1- 

^ 1 




JL^JLl'i' 

^ # ^^m^ ■ * ■ ■ 




, , , 1 , , , , 1 , , , , 1 , " 



0.5 1 1.5 2 2.5 3 

I A(t) 



Figure B.IO: Measured energy flow using calorimeter towers as a function of | A0| 
with respect to the jet direction. The measurements are compared to minimum 
bias Monte Carlo simulations. 

Energy Flow in dijet events 

A subsample of events is selected with at least two jets with pr > 7 GeV and \y\ < 
2.6, and the energy flow < >jcts are measured as a function of the rapidity 
separation between the two jets Ay^^ . For large Ay^^ , the measured energy flow, 
as defined above in the azimuthal direction, approximately decouples from the 
presence of a second jet and gives a more direct access to remaining underlying 
event contributions in the final state. Figure IB. 141 shows the measured energy 
flow compared to Monte Carlo simulations. As expected, for Ay^^ < 0.6 the 
presence of the second jet at |A0| ~ tt is very pronounced, while for Ay^^ > 1.2 
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Figure B.ll: Measured energy flow using calorimeter towers as a function of 
|A0| with respect to the jet direction, in different jet rapidity regions. The 
measurements are compared to minimum bias Monte Carlo simulations. 

a plateau of hadronic activity at large A(f) is clearly observed. Similarly, Fig. IB. 151 
presents the measurements based on tracks and only considering jets with \y\ < 
1.9. At large Ay^^ and large |A0| (in the plateau region) the different pythia 
simulated samples describe the data while the phojet sample underestimates 
the hadronic activity. 
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Figure B.12: Measured energy flow using tracks as a function of | A0| with respect 
to the jet direction. The measurements are compared to minimum bias Monte 
Carlo simulations. 
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Figure B.13: Measured energy flow using tracks as a function of | A0| with respect 
to the jet direction, in different jet rapidity regions. The measurements are 
compared to minimum bias Monte Carlo simulations. 
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B.2.1 Energy flow in the azimuthal direction 
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Figure B.14: Measured energy flow using calorimeter topotowers in dijet events 
and jets with pr > 7 GeV and \y\ < 2.6, as a function of | A0| with respect to the 
jet direction and the rapidity separation between the two jets. The measurements 
are compared to minimum bias Monte Carlo simulations. 
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Figure B.15: Measured energy flow using calorimeter topotowers in dijet events 
and jets with pr > 7 GeV and \y\ < 1.9, as a function of | A0| with respect to the 
jet direction and the rapidity separation between the two jets. The measurements 
are compared to minimum bias Monte Carlo simulations. 
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B.2.2 Energy Flow in rapidity 

The flow of energy around tlie jet direction is studied as a function of the distance 
in rapidity Ay to the jet axis. In this case only towers within a azimuthal 
distance A0 < D are considered (see Fig. lB.16|) . and energy instead of transverse 
momentum is computed to appreciate better the color flow between the jet and 
the direction of the proton remnants: 



< 



dE 
d(j)dy 



> 



jets- 



2L)|A?/|N 



^E(y-A|//2,0 + A|//2), 



(B.2) 



jet jets 



where E{y — Ay/2,(j) + Ay/2) is the sum of the energy of the constituents 
(calorimeter towers) at a given distance y to the jet. Bins of Ay = 0.2 have 
been considered. Measurements are performed for jets in different rapidity re- 
gions in the range —1.2 < y < 1.2 (see Figs. IB. 171 and IB. 181) and compared to 
Monte Carlo simulations. Assuming a symmetric response of the calorimeter for 
positive and negative rapidities, one can fold both figures in a single one and de- 
fine the observable in terms of absolute jet rapidity bins (see Fig |B.19| to reduce 
the statistical fluctuations. 
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Figure B.16: Sketch of the particle flow as a function of the distance in rapidity 
to the jet axis. 

The measured energy flow present the expected peak around the jet axis, and 
outside the jet cone the energy flow increases as the distance to the jet increases. 
The measurements show a structure in Ay that can be correlated with the transi- 
tion across different calorimeter subsystems with different response. The Monte 
Carlo simulation (PYTHIA MC09) provides a reasonable description of the core 
of the jet while systematically underestimates the activity away from the jet. 
The differences are particularly pronounced in the very forward region of the 
calorimeter and in the crack regions between calorimeters, and can be attributed 
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Figure B.17: Measured energy flow using calorimeter topotowers for jets with 
Pt > 7 GeV as a function of Ay in different jet rapidity regions. The measure- 
ments are compared to minimum bias Monte Carlo simulations. 

to deficiencies in the simulation of the detector material and response to low 
energy particles. 
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Figure B.18: Measured energy flow using calorimeter topotowers for jets with 
Pt > 7 GeV as a function of Ay in different jet rapidity regions. The measure- 
ments are compared to minimum bias Monte Carlo simulations. 



155 



Appendix B. Jet shapes and energy flow in pp collisions at ^/s = 900 GeV 



Beam 2009 



• 900 GeV MinBias DATA 
— 900 GeV Pythia MC09 



Anti-Kp D = 0.6, UNCORRECTED 
Statistical errors only 



o 



10 



i2 1 



A 
< 

Q 

(N 

i 

V 



10 



pi;' > 7 GeV 

; 1 y'^' 1 < 0.3 

ft ■ 

* ; 


: J : 

pf > 7 GeV 
0.3 < 1 y'" 1 < 0.6 


- 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 U 1 1 1 1 1 1 1 ■ 

pi;' > 7 GeV 

^ 0.6 < 1 y'"' 1 < 0.9 

: ^ ♦ : 

1 ... 1 ... 1 ... 1 ... 1 ... 1 ... 1 


p^^' > 7 GeV 
^ 0.9<|y'"|<1.2 

T ... 1 ... 1 ... 1 ... 1 ... 1 ... 1 



1 - 



6 -6 



Ay 



Figure B.19: Measured energy flow using calorimeter topotowers for jets with 
Pt > 7 GeV as a function of Ay in different jet rapidity regions. The measure- 
ments are compared to minimum bias Monte Carlo simulations. 
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B.2.3 Energy Profiles beyond the cone of the jet 

Finally, the energy profile < dp'r/dr'^ >iets is defined in a similar way as the 
differential jet shape but normalized to the area AA in each annulus, and without 
dividing, jet-by-jet, by the sum of the transverse momenta in < r < D: 

<jl^>^-ts=2x^T.P^i^-^^/'^^^ + ^^/'^)^ 0<r<D (B.3) 

Only central jets with rapidity in the region \y\ < 0.3 and \y\ < 0.6 are considered 
and the measurements are carried out using tracks. The use of tracks instead 
of calorimeter towers allows studying the particle flow around the jet direction 
without any bias due to the jet search algorithm employed, or the variations of 
the calorimeter response and the presence of calorimeter cracks as a function 
of rapidity. The measurements are compared to Monte Carlo simulations in 
Fig. IB. 211 As already mentioned, the Monte Carlo models overestimates the 
amount of transverse momentum close to the core of the jet. The different 
PYTHIA samples provide a reasonable description of the activity in the region 
away from the jet axis while PHOJET underestimates it. 




rapidity 



Figure B.20: Sketch of the particle flow as a function of the distance in radius 
to the jet axis. 
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Figure B.21: Measured energy flow using tracks for jets witli px > 7 GeV as 
a function of Ar for jets with \y\ < 0.3 and \y\ < 0.6. The measurements are 
compared to minimum bias Monte Carlo simulations. 
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The study of the energy flow away from the jet direction as described in Ap- 
pendix [B] was also performed in pp collisions at = 7 TeV at detector level. 

Figures fC . 1 1 and [C . 2 1 show the measured energy flow using tracks and calorime- 
ter clusters respectively, for jets with \y\ < 1.9 and 30 GeV < pr < 210 GeV. 
Figure [Cl3] extend in rapidity the measurement using clusters up to |y| = 2.8. As 
expected, the measured energy flow increases due to the presence of a second jet 
as A0 tends to vr. The measurements are compared to different MC predictions. 
PYTHIA-PERUGIA20101 gives a good description of the data. Both PYTHIA- 
MC09 and PYTHIA-DW tends to underestimate the activity in the region away 
from the jet axis, dominated by soft hadronic activity, whereas HERWIG++ 
tend to overestimate it, specially in the forward region. 
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Figure C.l: Energy flow using tracks as a function |A0| witli respect to tfie jet 
direction for jets with \y\ < 1.9 and 30 GeV < pr < 210 GeV. 
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Figure C.2: Energy flow using calorimeter clusters as a function |A(/)| with respect to 
the jet direction for jets with \y\ < 1.9 and 30 GeV < px < 210 GeV. 
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Figure C.3: Energy flow using calorimeter clusters as a function |A(/>| with respect to 
the jet direction for jets with 1.9 < \y\ < 2.8 and 30 GeV < pr < 210 GeV. 
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